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Abstract 

Neutrino oscillation data strongly support /i— r symmetry as a good approximate flavor symmetry 
of the neutrino sector, which has to appear in any viable theory for neutrino mass-generation. The 
/.J— r breaking is not only small, but also the source of Dirac CP-violation. We conjecture that both 
discrete /i— r and CP symmetries are fundamental symmetries of the seesaw Lagrangian (respected 
by interaction terms), and they arc only softly broken, arising from a common origin via a unique 
dimension-3 Majorana mass-term of the heavy right-handed neutrinos. From this conceptually 
attractive and simple construction, we can predict the soft fi^r breaking at low energies, leading 
to quantitative correlations between the apparently two small deviations ^23^ 45° and ^^13— 0° . 
This nontrivially connects the on-going measurements of mixing angle 623 with the upcoming 
experimental probes of ^13 . We find that any deviation of ^23 ~ 45° must put a lower limit 
on 013 . Furthermore, we deduce the low energy Dirac and Majorana CP violations from a 
common soft-breaking phase associated with fi—r breaking in the neutrino seesaw. Finally, from 
the soft CP breaking in neutrino seesaw we derive the cosmological CP violation for the baryon 
asymmetry via leptogenesis. We fully reconstruct the leptogenesis CP-asymmetry from the low 
energy Dirac CP phase and establish a direct link between the cosmological CP- violation and the 
low energy Jarlskog invariant. We predict new lower and upper bounds on the 613 mixing angle, 
1° < ^ 6° . In addition, we reveal a new hidden symmetry that dictates the solar mixing 
angle 612 by its group-parameter, and includes the conventional tri-bimaximal mixing as a special 
case, allowing deviations from it. 
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1. Introduction 



Neutrino physics has come to a golden era for a decade. The discovery of oscillations of solar and 
atmospherical neutrinos, and the further confirmations by using terrestrial neutrino beams produced 
in the reactor and accelerator experiments, have established the massiveness of neutrinos as well as 
the large mixing angles in the leptonic mixing matrix [ij^ . The two mass-squared differences (As, A^) 
and two large mixing angles (^12, ^23) are measured with good accuracy. A summary of the 3iy global 
analysis on these parameters [3] is shown in Table-lH with the conventions of A^ = m2 — rnl and 
Aa = \ml- {ml + ml)/2\ for the fit d [5] . 



Parameters 


A, (10- 


5eV') 


A, (10- 




sin^ (912) 


sin' 023 (023) 


sill' (013) 


Best Fit 


7.67 


2.39 


.312 (34.0°) 


.466 (43.0°) 


.016 (7.3°) 


Icr Limits 


7.48- 


7.83 


2.31 - 


2.50 


.294- .331 
(32.8° - 35.1°) 


.408 - .539 
(39.7° -47.2°) 


.006 - .026 

(4.4° -9.3°) 


2cr Limits 


7.31 - 


8.01 


2.19- 


2.66 


.278- .352 
(31.8° - 36.4°) 


.366 - .602 
(37.2° -50.9°) 


< .036 
(< 10.9°) 


3cr Limits 


7.14- 


8.19 


2.06- 


2.81 


.263 - .375 
(30.9° - 37.8°) 


.331 - .644 
(35.1° - 53.4°) 


< .046 
(< 12.4°) 


90%C.L. Limits 


7.36- 


7.94 


2.24- 


2.60 


.283 - .344 
(32.2° - 35.9°) 


.380 - .582 
(38.1° -49.7°) 


< .032 
(< 10.4°) 


99%C.L. Limits 


7.23- 


8.13 


2.10- 


2.75 


.270- .365 
(31.3° - 37.2°) 


.344 - .627 
(35.9° - 52.3°) 


< .042 
(< 11.8°) 



Table 1: The global 3z^ fit [3] for the neutrino mass-squared differences and mixing angles including 
the available data from solar, atmospheric, reactor (KamLAND and Chooz) and accelerator (K2K 
and MINOS) experiments [5]. 

The Table-[T] shows that the neutrino mixings pose a striking pattern with two large mixing angles 
{9i2, ^23) and a small mixing angle ^13 , very different from that of quark sector. It is also intriguing 
to note that the central value of ^23 is already somewhat below the maximal mixing angle[^ (45°) and 
613 is favored to deviate from zero at 95%C.Li To be specific, the data allow two small deviations 
of the same order, at the 2a level, 

- 7.8° < (023 - 45°) < 5.9° , 0° ^ (013 - 0°) < 10.9° , (1.1) 

with the best fitted values, (023 — 45°) = —2.0° and (0i3 — 0°) = 7.3°. This naturally provides a 

^The leptonic mixing matrix is conventionally called Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [2] whose 
form will be discussed later in Sec. 13. II 

^It was noted back in 2004 [7] that the data have slightly favored a smaller 623 than its maximal value. Although this 
intriguing small deviation 6^23 — 45° < is not yet statistically significant, the tendency of having a negative 623 — 45° 
is found to be robust, due to the excess of e-like atmospheric neutrino events in the sub-GeV sample [7]. 

very recent oscillation analysis on ^13 gives, sin^ ^13 = 0.02 ± 0.01 at la level, with the central value slightly 
shifted upward [5] . 



3 



fairly good zeroth order approximation, O23 = 45° and ^13 = 0°, under which two exact discrete 
symmetries emerge, namely, the /i— r symmetry [TTj and the Dirac CP conservatior0 in the neutrino 
sector. We stress that the /U — r symmetry (as well as associated Dirac CP-invariance), as a good 
zeroth order approximation reflected by all neutrino data, has to appear in any viable theory for 
neutrino mass-generation. 

Since the /x— r symmetry, as a good zeroth order flavor symmetry of the neutrino sector, leads to 
^13 = and thus the Dirac CP-conservation, the /i — r symmetry breaking terms must be small and 
also serve as the source of the Dirac CP-violation. On the theory ground, it is natural and tempting 
to expect a common origin for all CP-violations, although the Dirac and Majorana CP-violations 
appear differently in the low energy effective theory of light neutrino mass-matrix. Given such a 
common origin for two kinds of CP-violations, then they must vanish together in the h—t symmetric 
limit. 

With these key observations, we conjecture that both discrete /U — r and CP symmetries are 
fundamental symmetries of the seesaw Lagrangian (respected by interaction terms), and they are 
only softly broken, arising from a common origin via a unique dimension-3 Majorana mass-term of 
the heavy right-handed singlet neutrinos. The reason for the fJ, — T and CP breaking terms being 
soft is because we consider both symmetries to be respected by interaction terms of the seesaw 
Lagrangian |9j , thus the unique place for such breakings is the dimension-3 singlet Majorana mass 
term of right-handed neutrinos. 

From the above conceptually attractive and simple construction with the soft breakings of two 
discrete symmetries from a common origin, we can predict the soft fJ, — T breaking at low energies, 
leading to quantitative correlations between the apparently two small deviations U, = ^23 ~ 

45° 

and 5x = ^13— 0°. The mixing angle ^13 will be more precisely probed at the upcoming reactor 
experiments Double Chooz[lJ], Daya Bay^], RENO [16] as well as the accelerator experiments 
T2K [18] and NOi^A [T7], etc|j, while improved measurements of the mixing angle 623 are expected 
from the Minos [20] and T2K [18j experiments etc. The future neutrino factory and super-beam 
facility [2T] will further pin down these key parameters with high precision. Finally, we further 
derive the cosmological baryon asymmetry rj^ via leptogenesis and analyze the interplay between the 
leptogenesis scale Mi and the low energy Jarlskog invariant J as well as neutrinoless double-beta 
decay observable M^e- 

Table-|T] also shows that the solar angle 9i2 is best measured among the three mixing anglej^, 



^Here the possible Majorana CP-phases may still appear, but in our theory construction fSec. l2.2 |l. they are not 
independent and will vanish as the Dirac CP violation goes to zero in the /i — r symmetric limit. 

^It is interesting to note that a nonzero 023~45° is further motivated by the quark-lepton complementarity [TD] via the 
relation + 623 = 45°, where experiments have measured the quark-sector CKM (Cabibbo-Kobayashi-Maskawa [12]) 
mixing angle ~ (2.36 ±0.06)° [13], causing a deviation in the leptonic mixing angle 623 — 45° — 6'23 < 45°, consistent 
with the data in Table-[T] 

®For instance, the proposed LENA experiment [19] will also probe 613 and CP-violations with good precision via 
long baseline neutrino oscillations from CERN to the LENA detector at Pyhasalmi mine (2288 km apart). 

^Table-[T] shows that the so-called tri-bimaximal (TBM) mixing ansatz, sin^ 6^12 = | (S12 = 35.3°), is above the la 
upper bound on 6^12. A latest three-flavor oscillation analysis of SNO Collaboration [22] combined data from all solar 
experiments and KamLAND reactor antineutrino experiment [23] . this gave an even tighter constraint on the mixing 



4 



so it is not our major concern. The main goal here is to examine the least known angle ^13 (to 
be measured at Double-Chooz [H], Daya Bay [E], RENO [BJ, T2K [Igj and NOi/A [IT] experiments) 
which is crucial for nonzero Dirac CP-violation, as well as its quantitative correlation with the 
deviation of ^23 — 45° from the naive maximal mixing (which will be further probed by MINOS |20j 
and T2K [T8] experiments, etc). We can predict nontrivial lower and upper limits on the mixing 
angle 6*13. 

This paper is organized as follows. In Sec. [2] we present a unique construction of the soft fi—r and 
CP breakings from a common origin in the neutrino seesaw. Then, in Sec. [3] we present a low energy 
reconstruction of the light neutrino mass matrix with /i — r and CP violations including the small 
parameters 5a and 6x ■ With these, we will further derive, in Sec.|H the low energy /i — r and CP 
violation observables from the common soft breaking in the neutrino seesaw. As the 6a and 6x arise 
from the common origin, they are proportional to a (small) common soft breaking parameter at the 
next-to-leading order of the well-defined expansion. So we can make quantitative predictions for their 
correlations, and discuss their implications for the upcoming experimental probes. Then, in Sec.[5]we 
analyze the cosmological CP violation via leptogenesis |24|, [25] in our model, which provides the origin 
of matter — the baryon asymmetry of the universe. After inputting all the known neutrino data 
and the observed baryon asymmetry [26|. [13] , we can derive the direct link between the cosmological 
CP- violation and the low energy Jarlskog invariant [27] . We further place a lower bound on the 
leptogenesis scale for producing the observed baryon asymmetry, and deduce new lower and upper 
limits on the mixing angle ^13 . In Sec.[6l we generally prove why the solar mixing angle does not 
receive any correction from the soft fi—T and CP breakings. We show that 9 12 is actually protected 
by a new hidden symmetry in the seesaw Lagrangian, and its value is dictated by the group parameter 
of this hidden symmetry. We finally conclude in Sec. [71 

2. Soft ^i — T and CP Breaking Originated from Neutrino Seesaw 

We conjecture that both the discrete /i — t and CP symmetries are fundamental symmetries of 
the neutrino seesaw Lagrangian (respected by interaction terms), and are only softly broken from a 
common origin. In fact, the only place in the seesaw Lagrangian which can provide such a common 
origin is the dimension-3 singlet Majorana mass-term of right-handed neutrinos. The Dirac mass- 
term cannot provide such a soft breaking as it is generated by the dimension-4 Yukawa interactions 
after spontaneous electroweak symmetry breaking. In Sec. l2.1[ we first consider the minimal neutrino 
seesaw Lagrangian with exact /i — r and CP invariance, and derive the seesaw mass-matrix for the 
light neutrinos (from which we deduce the zeroth-order mass-eigenvalues and mixing angles). This 
/i — r and CP symmetric limit predicts the mixing angles (^23) ^13) = (45°, 0°). Then, in Sec. [221 
we will construct a unique soft breaking term providing a common origin for both /U — r and CP 
breakings. Prom this we will further derive predictions for the /U — r and CP breakings in the low 
angle, 612 ~ 34.38lQg7 (degrees), at la level, which also leaves the TBM mixing value around the la upper boundary. 
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energy light neutrino mass-matrix, by treating the small soft-breaking as perturbation up to the first 
nontrivial order (SecU]). The predictions for the seesaw-scale leptogenesis and its correlations with 
low energy observables will be analyzed in SecO 



2.1. Neutrino Seesaw with /x — r and CP Symmetries 

Seesaw mechanism [S] provides a natural explanation of the small Majorana masses for light neutrinos 
which violate lepton number by two unitqj. For simplicity of demonstration, we consider the La- 
grangian for the minimal neutrino seesaw [291 130j . with two right-handed singlet Majorana neutrinos 
besides the standard model (SM) particle content, 

Css = -LYiMji - LY^'^M+^J^MRCAf + h.c. 

= -II MiiR - ulmj^N + ^A/^MrCM + h.c. + (interactions) , (2.2) 

where L denotes three left-handed neutrino-lepton weak doublets, i = (e, fi, t)'^ contains charged 
leptons, Uj^ = (z^g, u^, v^y is for the light flavor neutrinos, and M = (A'^^, A',-)-'" represents two 
heavy right-handed singlet neutrinos. We have also denoted the SM Higgs doublet by $ {^) with 
hypercharge \ {—\)-, and the 3x3 lepton and neutrino Yukawa-coupling matrices by Y^ and Y^, 
respectively. The lepton Dirac-mass-matrix Mi = vYil\f2 and the neutrino Dirac-mass- matrix 
vTLj^ = vY^j^pl arise from the dimension-4 Yukawa interactions after the spontaneous electroweak 
symmetry breaking, (<1>) = (0, )"^ 7^ , and the dimension-3 Majorana mass-term for Mr is a 
gauge-singlet. 

This minimal seesaw Lagrangian in Eq. ()2.2p can be regarded as an effective theory of the general 
three neutrino seesaw where the right-handed singlet is much heavier than the other two (A^^^j, A'^,-) 
and thus can be integrated out at the mass-scales of (Np,, N^), giving rise to Eq. (j2.2|) . Consequently, 
the minimal seesaw generically predicts a massless light neutrino p9] , which is always a good approx- 
imation as long as the lightest left-handed neutrino has its mass much smaller than the other two 
(even if not exactly massless). As we will prove in Sec. 6.3, the general three-neutrino-seesaw under 
the /u— T symmetry (^2^) and hidden symmetry (Z|) also predicts a massless light neutrino, sharing 
the same feature as the minimal seesaw we consider here. 

After integrating out the heavy neutrinos (A'^, A^^^) from Eq. ()2.2p . we can derive the seesaw 
formula for the 3x3 symmetric Majorana mass-matrix of the light neutrinos, 

M^ ~ m£,M^V|^. (2.3) 

In general it is expected that the lepton sector would obey a flavor symmetry Ge different from the 
H — T symmetry in the neutrino sectoip, due to the large mass-splitting between (/i, r) leptons. 
The two symmetries and Gi could result from spontaneous breaking of a larger flavor symmetry 



*It was shown [JT] that the seesaw mechanism can also be reahzed for generating small Dirac masses for light 
neutrinos with lepton-number conservation. This direction will not be explored in the current study. 
^The A* — T symmetry Zj^ will be defined shortly, in Eqs. H2.5[l - ()2.6p and the text just above them. 
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Gp [28]. The invariance of lepton mass matrix under the transformation of left-handed leptons 
(G G^) is gI^M^m}g^^ = M^Mj . The mass matrix MiMj can be ensured to be diagonal via 



proper model-buildings, e.g., via making diagona 
we can make the simplest choice of = Z3 , witlj^ 



and nondegenerate. For an explicit realization, 



GiL = diag(l, rj, rf) , (2.4) 

where rj = e^'^'^l'^ and G\^ = I with / being unit matrix. Thus, the above invariance equation 
G^j^M^M^G^^ = M^Adj makes M^M^ fully diagonal, and corresponds to the left-handed leptons in 
their mass-eigenbasis. This means that the PMNS matrix V = ujj^U^^ = U^^ . The right-handed 
leptons can be further rotated into their mass-eigenbasis, without affecting the PMNS matrix, except 
having a diagonal lepton-mass-matrix Mi in Eq. (j2.2p . (Note that the relation V = U^^ holds as 
long as MiAlj is ensured to be diagonal via the above symmetry transformation Gil £ , but this 
does not require the lepton mass-matrix Mi itself being generally diagonal because the right-handed 
leptons can still be in the non-mass-eigenbasis.) 

The Lagrangian (|2.2p is defined to respect both the — r and CP symmetries. So the Dirac 
and Majorana mass matrices m^-, and are real, as well as symmetric under the transformations, 
i/^ -H- pUr and -H- p'Nr , where p,p' = it denote the even/odd parity assignments of the light 
and heavy neutrinos under the /U — r symmetry Z2 ^ . This means that the mass matrices and 
Mj^ obey the following invariance equations, 

flm^fji = mj,, T^Mj^Tr = Mji, (2.5) 

with 

where we have added a "tilde" on and M/j to indicate their most general forms including possible 
T parities, and later we will remove the "hat" after some useful simplifications. Thus we can deduce 
and to have the following structure, 

/ a p'a \ 
fflD = be, 
\ pp'c pp'b J 

which are all real due to the CP conservation. Here, the hatted quantities denote the mass matrices 
before rephasing, while M22 and M23 denote the absolute values of the {22} and {33} elements of 
with their signs defined as and P/jP/j, respectively. Note that the Dirac mass-matrix contains 
only three independent real mass-parameters (a, b, c), and the symmetric Majorana mass-matrix Mr 




Mr=Pr[,^^ 1 |, {p,p\pr,p'r = ±), (2.7) 
VrM2 



'23 



M22 



^''Actually any discrete group Z„ with n J5 3 can provide a diagonal nondegenerate 3-dimensionaI representation 
and thus does the same job [2H] • 
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has just two, {M22, ^^23)- Substituting the above and into the seesaw equation ()2.3p . we 
deduce the mass matrix for hght neutrinos, 



2a^ 



PR 



A/22 +p'p^ A/23 



v 



a{b+p'c) 
M22+p'p^A/23 



A/22 A/23 



M22-p'p^A/23 



M22+p'p^Af23 A/22-p'p^A/23 



pa(b+p'c) 
A/22+p'p^A/23 

(fe+p'c)^ (b-p'c)^ 



A/22 +pVfl A/23 A/22 -p'pij A/23 



which we call the zeroth order mass-matrix since we will further include the small soft-breaking effect 
in the next subsection. 

For the mass formula (j2.8p . some comments are in order. First, it is clear that the sign factor 
p' = ziz (due to the parity of the /U— r transformations under Zg^) accompanies c and M23 everywhere 
in Eq. (|2.8p . Besides, is always associated with M23, as expected. It is convenient to simply rescale 
c as p^c; then p' and p'j^ always accompany each other in (j2.8p , and appear together in front of both c 
and M23. Second, we note that the {12} and {13} elements of Mi, only differ by an overall sign p = it 
(due to the other parity factor under Zg^); the sign p also determines the sign of the mixing angle 
^23- As our sign-convention, we always define 623 £ [0, ^] ^ 0, so we can make a simple rephasing 
for the light neutrino u^, i.e., (u^, f^, u^)'^ — )• Vhi^e^ i where Vl is a 3 x 3 diagonal matrix, 

Vl = diag(l, l,p) . This will not affect the leptonic PMNS mixing matrix since it only contributes 
to the rephasing matrix U" as defined in Sec. 13. II (which is not related to the PMNS matrix). Under 
this rephasing, the mass matrix ()2.8p becomes Vj^MyVi-]!. Third, we note that the overall sign pp^ 
in ()2.8p originates from the overall sign in (|2.7p and can be rotated away by uniformly rephasing the 
Majorana neutrino fields N — t- ^/P^J^ ', we can further remove the relative sign p^ in M/j [ ()2.7p ] by 
a rephasing of M^i or A/V- In summary, to simplify the sign-conventions we just need to make the 
following combined transformations for light and heavy neutrinos, in one step, 



with Vp = ^/^diag(l, 1, p) , = ^/;p^diag(l, p'^) , where v = (z/^, i/^, v^) 
mass matrices in (|2.7p transform accordingly. 



with 



m 



D 



Mr 




where 



P = p'p'r = ± 



R 



Mr 

'M22 M23' 
^M23 M22, 

M23 



Mr = V'rMrVr, 



(2.9) 

Then the two 
(2.10) 



M: 



22 



M 1/ 



(2.111 



l-r > 0. 



22 



(2.12) 



^'^From Eq. (|2.8|) or Eq. 1)2. 16p below, we have det(M^) = det(Af^) = 0, which actually holds in general minimal 
seesaw (even after including p — r and CP violations). 
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In the above we have simply rescaled c as p^c (which does not matter since we can fully eliminate 
the parameter c in our final relations between physical observables) . Since both M22 and M23 denote 
absolute values, we always have R = 1 — r > , i.e., r < 1 . 

Furthermore, the r symmetry transformation matrices Tl and Tr will change into Tl and Tr 
accordingly, 



with 



Tl 




Tr 




(2.13) 



(2.14) 



Thus, the fJ. — T invariance equations for the Dirac and Majorana mass-matrices in ()2.1ip become. 



Tlm^Tj^^ 



m 



D 1 



T^MrTj^ = Mr, 



(2.15) 



which can be readily verified. Finally, we rewrite the symmetric seesaw mass-matrix M^, in ()2.8p as. 



2a^ 



22 



l+pR 



a{b+pc) 
1+pR 

{h+pcf {b-pcf 
l+pR l~pR. 



\ 



a{b+pc) 
1+pR 

(b+pc)^ _ (b-pc)^ 
1+pR 

{b+pcf 
l+pR 



\ 



+ 



l-pR 

{b~pc)^ 



l~pR 



(2.16) 



which is invariant under /x — r symmetry, 



tIm^Tl 



(2.17) 



As mentioned earlier, we can readily verify det(Mj,) = 0, which actually holds for any two- neutrino 
seesaw. 

Diagonalizing the /x — r and CP symmetric mass-matrix (j2.16p . we derive the mass eigenvalues 
and mixing angles 

2a'^ + {h + pcf 



mi 



0, 



tan 9i2 



rrio 



V2\c 



M22\l+pR\ 



723 



45° , 



rrin 



{b - pcY 
M22\l-pR\ ' 



0°, 



(2.18a) 



(2.18b) 



\b + pc\ 

where we have made all mass-eigenvalues positive and the three mixing angles within the range 
[0, ^] by properly defining the rotation matrix (cf. Sec. [3] and Sec. l6.ip . Here the mixing angles 
(^23) ^13) = (45°, 0°) are direct consequence of the /x — r symmetry, while the value of 612 does not 
depend on it. The determination of O12 from additional flavor symmetry will be further analyzed 
in Sec. El The Eq. ()2.18ap shows that the mass-spectrum of light neutrinos falls into the "normal 
hierarchy" (NH) pattern (mi < m2 < m^). From (|2.1ip . we can diagonalize the /i — r and CP 
symmetric mass-matrix Mr for right-handed neutrinos, with two mass-eigenvalues given by 



Ml 



\r\M' 



22 



M2 = (2-r)M22, 



(2.19) 
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where 2 — r > is ensured because of R = 1 — r > and thus r < 1 . 

Finally, we comment on the sign choice for p = zt . From Table-dl we deduce, at 2a level. 



1712 / ^"^^21 I ^ 



^ = 0.17-0.19 < 1. (2.20) 



So we can derive the mass ratio, 

7712 20^ + {b + pcf \1 - pR\ 



<. 1, (2.21) 



7Tt3 {h — pcY \\ + pR\ 

where we note th,^ = as given in ()2.18ap (which holds even after including /i— r and CP violations). 
Using the definition i? = 1 — r > and the mass-eigenvalues in ()2.19p . we can rewrite the ratio, 

I\A Ml ,^ _ , 

= = = f? , (tor p = + , 

2-r M2 ^' V 
(2.22 
2-r M2 _i _ , 

= = Mi='^ ' (forp = -). 

The Dirac mass-matrix raj^ in ()2.1ip arises from the products of Yukawa couplings and Higgs vacuum 
expectation value. So for a natural seesaw without fine-tuning where Yukawa couplings are of 0(1), 
we expect tti^ = 0(100 GeV) , lying at the weak scale, so the ratio (2a^ + (6+j3c)^)/(6— pc)^ = 0(1) . 
Hence, we deduce from ([2:21]) and ([222]), 

Ml 

g = — i « 1 , (for p = +) , (2.23a) 

M2 

= ^ « 1, (forp = -). (2.23b) 

This indicates that the mass-difference between the heavy neutrinos is large, which is also what we 
need for the leptogenesis analysis in SecO From (|2.22p we can further resolve the parameter r , 

= 2p,g + 0(g2) « 1, (forp = +), 

1 + pM 

(2.24) 

= 2(l-p,(?-i) + 0(g-2) ~ 2, (forp=-), 

where = ± denotes the sign of r. We see that the case p = — leads to r ~ 2 which contradicts 
with the original condition i? = 1 — r > . Hence, the sign-choice of p = — is excluded (or strongly 
disfavored), and we can focus on the case with p = + from now on. 



2.2. Common Origin of Soft fJ, — T and CP Breaking in Neutrino Seesaw 
2.2.1. Unique Construction of Soft /x — r and CP Breaking in Neutrino Seesaw 

As we pointed out in the introduction (Sec.[T]), the ^—t symmetry serves as a good zeroth order flavor 
symmetry of the neutrino sector, it further leads to ^13 = and thus the Dirac CP-conservation. 
This means that the fi — r symmetry breaking is not only small, but also the source of the Dirac 
CP-violation. On the theory ground, it is natural and tempting to expect a common origin for all 
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CP- violations, although the Dirac and Majorana CP-violations appear differently in the low energy 
effective theory of light neutrino mass-matrix. With such a common origin for two kinds of CP- 
violations, then they must vanish together in the /i — r symmetric limit. 

With the above key observations, we conjecture that both discrete /i — r and CP symmetries are 
fundamental symmetries of the seesaw Lagrangian (respected by interaction terms), and are only 
softly broken from a common origin in the seesaw sector. We observe that the only place for such 
a soft breaking is the dimension-3 singlet Majorana mass-term of the right-handed heavy neutrinos 
M in Eq. (j2.2p . The Dirac mass term cannot provide such a soft breaking as it is generated by the 
dimension-4 Yukawa interactions after spontaneous electroweak symmetry breaking. 

Under this conjecture, we can thus define the common soft-breaking term for ^ — t and CP 
symmetries in the Majorana mass- matrix Mr , via a single complex parameter ^e*'^ , in a unique 



where the module < C < 1 and the phase angle £ [0,27r). Since the sign of C, can always be 
absorbed into the phase e*"^, we only need to consider > . Besides, C, should be significantly 
smaller than one as the /i — r and CP violations are small. As explained in Sec. ()2.ip . the elements 
M22 and M23 will be taken in the positive range without losing generality. We note that the unique 
CP phase e*'^ in (|2.25|) will not only generate the low energy CP-violations (via both Dirac and 
Majorana phases) appearing in the neutrino oscillations and neutrinoless double-/? decays, but also 
provide the CP asymmetry needed for the successful leptogenesis of baryon asymmetry (cf. Sec lS.ip . 

We also note that this soft-breaking term cannot be located in the {23} and {32} elements of 
Mn because the Majorana mass-matrix Mr is symmetric and the equality between the two off- 
diagonal elements automatically respects the /x — r symmetry. Hence, the only place to softly and 
simultaneously break both /ir and CP symmetries is in the diagonal elements of Mr. Then, we may 
wonder whether we could relocate the soft-breaking term in the {ll}-element instead of {22}-element 



But we can readily prove that the definition ()2.26p is actually equivalent to that in Eq. ()2.25p for our 
present study (where we are only interested in deriving correlations among the physical observables). 



way. 




(2.25) 



of Mr 




(2.26) 



To see this, we can simply rename the right-handed neutrinos N 



{N^, N^) as Af' = {N^, , 



I.e., 




with 




(2.27) 
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Then, the mass matrices rrij^ in ()2.1ip and Mr in ()2.26p change accordingly, 

/pa a\ /a' pa'\ 
nijy m'^ = mnTZ = c 6 = 6' c' , (2.28a) 



/pa 


a) 




fa' 


pa'\ 




6 




b' 




\pb 


pcj 






pb'l 



, ^ /M22 M23 \ 

Mr Mr = TZ^MrTI = , (2.28b) 

^ \M2, M22{l-Cen) 

where we have simply renamed, (pa, c, 6) = (a', 6', c') . From this we see that the Majorana mass- 
matrix m'r is identical to Mr as defined in Eq. (|2.26p , while the Dirac mass- matrix takes exactly 
the same form as in Eq. ()2.1ip except all its elements are primed in the new notation. Since our 
final physical results in Sec. 1310 only concern the correlations among physical observables and do not 
explicitly rely on either (a, 6, c) or (a', 6', c')0, we conclude that the parametrization in (j2.26p is 
physically equivalent to ()2.25p for our study, thus the formulation of the common origin for soft /x— r 
and CP violations is indeed unique under our construction. In consequence, we only need to focus 
on the analysis of Eq. (|2.25p for the rest of our paper. 

In summary, the Dirac mass- matrix and Majorana mass-matrix Mr (with soft /U — r and CP 
breakings from a common origin) can be uniquely parameterized as follows, 

I " \ /M22 M23\ 

= 6 c , Mr =[ , M33 = M22 (1 - Ce*") , (2-29) 

t) V^23 M33y 

where we have set the sign p = + for as discussed by the end of Sec. l2.1[ 
Thus, we can explicitly derive the seesaw mass-matrix for light neutrinos. 

My = niDM^^rnj^ 

/a2(M22 + M33-2M23) a[cM22 + 6M33-(6-hc)M23] a[6M22-hcM33-(6-hc)M23] \ 

1 

c2M22 + 62M33-26cM23 6c(M22 + M33) - (6^ + c2)M23 

\ 62M22 + c2M33-26cM23 / 



IIM^II 




(2.30) 

where we have denoted, \\Mr\\ = det(Mij) = M22M33 - Mfg . We find that it is useful to further 
decompose Mi, into /i— r symmetric and anti-symmetric parts. 

My = M^ + M^, (2.31) 
/A Bs Bs\ 

Mt = - (M, + tIm.Tl) =1 Cs dV (2.32a) 

Csl 



^^We note that under our present construction these parameters wiU aU be ehminated from our final results since 
they are neither direct observable nor theoretically predicted (cf. Sec. l4l5p . 
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with 



and 



M« = ^ (M, - tIm.Tl) =1 Ca I , (2.32b) 



M22+M33-2M23 2 .0 QQ ^ 

^ = " ' (2.33a) 

^1+^2 M22 + M33-2M23 

= = ^py^^^ «(^ + -)' (2-33b) 

r^^n -1 ' , ^ M22+M33-2M23 , ,2 .oQQ^ 

C. + i^ = = ^p^^ {b + c) , (2.33c) 

2\\Mr\\ V ; , V ; 



n — XT3 B1-B2 M22-M33 

^ — xr' C*! - M22-M33 2 ,2n Q/|U\ 

^» = '^^'^ = = 2||Mh|| ('-^'^^ 





2 




+ C2 




2 




+ C2 


2 



where 



= Ce*" = 0(C), (2.34c) 

M22 

showing that the parameters 5Ba and 5Ca must vanish in the /i— r (and CP) symmetric hmit ^ — )• . 
So the ^i — T anti-symmetric part is a net measure of the /i — r breaking. From (|2.34p we further 
compute the ratio of 5Ba and 5Ca, 

TTT = T"^ — ' (2-35) 

OCa + C 

which is a real number, independent of the /i— r breaking and depending only on the elements of the 
/i— r symmetric rrij^ . Note that in (|2.34|) . all fi—T breaking effects are contained in the coefficient, 
(M22 — M33)/||M/j|| , which is exactly canceled in the ratio 6Ba/6Ca- In Sec. l6.Hl6T2l we will further 
prove that the ratio (|2.35p is directly connected to tan^^ for the solar mixing angle. 

2.2.2. Soft fi — T and CP Breaking via Perturbative Expansion 

According to our construction in Sec. 12. 2. 11 the size of the soft ^u— r and CP breakings are characterized 
by the small quantity \(\ < 1 , which can be treated as a perturbative expansion . As we noted by 
the end of Sec. 12. 11 there is another small /U — r and CP symmetric parameter |r| < 1 , appearing in 
the Majorana mass-matrix of heavy right-handed neutrinos, 

/I 1-r \ / i 

Mr = M22 = Mio ( ^ , (2.36) 
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where we have used the notations, 



r = I- R = 1 



M23 
M22 



Al- 



io 



rM22 , X 



Oil). 



(2.37) 



With these we can exphcitly write the hght neutrino mass matrix (j2.30p as, 



/ra?{2-X) ra{b+c-bX) 



ra 



(b+c-cX) \ 



{2-r-X)M: 
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V 



(2.38a) 



{b-cf+rb{2c-bX) -{b-cf +r{b'^ -bcX) 

{b-cf+rc{2b-cX) J 

Since the unique /u — r and CP violation quantity Qe^^ is now fully hidden in the ratio X = ^e*"^ = 
0(1) , we see that the only superficial small parameter for the effective perturbative expansion is 
r . So, we can reorganize the light neutrino mass-matrix in terms of the r-expansion alone, but 



with X unexpanded, 



where 



{b-cf 



{2-X)A'ho 



(2-X)2Mio 





1 -1 
1 

^2„2 



/{2-Xya' {2-X)[{l-X)b+c]a {2-X)[b+{l-X)c]a\ 
[(l-X)6+c]2 {l-X){b+c)'^+X%c 

[b+{l-X)c]' J 



\ 



(2.39) 



(2.40a) 



(2.40b) 



We note that the zeroth order mass-matrix Mj^*^^ is complex, but the CP-phase only appears in 
an overall factor and thus causes no observable CP-violation. The above decomposition can be 
symbolically denoted as. 



A Bi B2 
Ci D 
C2 



^ Bq Bq 

Co Do \ + 
\ Co 



/ 5A 6B1 6B2 \ 

6C1 6D 
\ SC2 J 



(2.41) 



with its zeroth order matrix corresponding to ()2.40ap and the next-to-leading order (NLO) matrix 
equal to (|2.40bp at 0{r). 

For convenience, we further decompose 5Mi, into the fi — T symmetric and anti-symmetric parts, 

/ 5Ba -5Ba 



6M^P = 6M' + 6M, 



6A 5Bs 5Bs \ 
6Cs 5D 

SCs J 



+ 



V 



6Ca 

-SCa 



(2.42) 



with 



5B, = 



5Bi + 5B2 



SBa = 



6B1 - 6B2 



(2.43a) 
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6C, = 



5Ci — 6C2 



(2.43b) 



2 ' 2 

This decomposition is actually unique. For our current model with the expansion up to 0{r), we 
deduce from ^^IM^ - ^im and (IXi2D - (ITi5D . 



An 



5n 



Cn 



and 



SA 

5Bn 



Mio ' 
a{b + c) 
2Mio ' 
a(6 - c)X 
~ 2(2-X)Mio 



51? 



^° (2-X)Mio ■ 

[(fe+c)^(l-X) + bcX^] 

(2-X)2Mio 
[2(6+c)^(l-X) + (6^+c^)X^] 
2(2-X)2Mio 
(62-c2)X 



2(2-X)Mj 



r . 
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We find it useful to further decompose the /i — r symmetric part into, 



/6A 5Bs 5Bs 

6Cs 6D I = 



(5A 5Bs 6Bs 

SEs 6Es I + 



/O 



6E' -5E\ 




with 



8E, 



5E' 



- {5Cs + 6D) 



-{5Cs-SD) 



{b+cf 
4Mio 



4(2-X)2Mio 



(2.44) 



(2.45) 



(2.46) 



(2.47a) 
(2.47b) 



Finally, we note that from ()2.45p we can compute the ratio, 

SBa a 



(2.48) 



5Ca b + c' 

which only depends on the elements of the /i — r symmetric Dirac mass-matrix m^-,, and agrees to 
Eq. (j2.35p at the end of Sec. 12. 2.1] where this ratio was generally derived without invoking the explicit 
form of Adji and without making the r-expansion. As we will show in Sec. l6.Hf6.21 this ratio is related 
to the solar mixing angle via tan^^ and is independent of the /i— r breaking from the Mr. 

Before concluding this section, for clarity we summarize the key features of our approach together 
all assumptions we have used. As explicitly stressed in Sec. 1 and at the beginning of Sec. 2, the 
oscillation data (Table-1) strongly support the /i — r symmetry as a good approximate symmetry of 
the neutrino sector, rather than the lepton sector (due to the large mass-splitting between the (/i, r) 
leptons). Thus the natural realization of /U— r symmetry is in the diagonal basis of left-handed leptons 
(corresponding to a diagonal mass- matrix product M£Mj ), which can be easily realized by a proper 
lepton flavor symmetry like Z3 [cf. Eq. ()2.4p ]. so the measured PMNS mixing matrix V = Uyi 
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arises from diagonalizing the light neutrino mass-matrix. Because the fi—T symmetric hmit enforces 
^13 = and thus vanishing Dirac CP-violation in the PMNS matrix, the Dirac CP-violation must 
originate from the fi—r breaking. On the theory ground, it is tempting to expect a common origin for 
all CP-violations (although the Dirac and Majorana CP-phases appear differently in the low energy 
light neutrino mass-matrix). With these key observations, we have conjectured that both discrete /i— r 
and CP symmetries are fundamental symmetries of the seesaw Lagrangian (respected by interaction 
terms), and they are only softly broken, arising from a common origin; the only place for such a 
soft breaking is the unique dimension-3 singlet Majorana mass-term for right-handed neutrinos. This 
means that Dirac mass-matrix naturally respects /x— r and CP symmetries, and only the singlet 
Majorana mass-matrix Mji can softly break them. This conjecture can hold for general neutrino 
seesaw with three right-handed neutrinos. In our explicit formulation, we have chosen the minimal 
neutrino seesaw |29t [50] (with two right-handed neutrinos added to the SM particle content) and 
proved that the common origin of the soft ^—t and CP breaking is uniquely characterized by a single 
complex parameter in M/j [cf. Eg. ()2.29p in Sec. 2.2.1]. As explained above Eg. (|2.3p . the minimal 
seesaw [291 l30j can be viewed as a useful effective theory of the general three neutrino seesaw with one 
right-handed neutrino much heavier than the other two and thus being integrated out of the seesaw 
Lagrangian (j2.2|) . It serves as a good approximation as long as the lightest left-handed neutrino has 
its mass much smaller than the other two (even if not exactly massless). We will further prove in 
Sec. 6.3 that the general three-neutrino-seesaw under the /x— t symmetry (Zg^) and hidden symmetry 
(Z|) also predicts a massless light neutrino, sharing the same feature as the minimal seesaw. 



3. Model-Independent Reconstruction of Light Neutrino Mass 
Matrix with (jL — t and CP Violations from Low Energies 

In this section, we present the model-independent reconstruction of the Majorana mass-matrix for light 
neutrinos, in terms of the low energy observables (mass-eigenvalues, mixings angles and CP phases). 
This reconstruction will be effectively expanded by experimentally well-justified small parameters up 
to the next-to-leading order. Then, in the next section we can readily apply this to our soft-breaking 
model (defined in Sec. l2.2p . and systematically derive our model-predictions for the correlations among 
the low energy observables as well as for the link to the leptogensis at seesaw scale. 



3.1. Model-Independent Reconstruction Formalism 

The 3x3 Majorana mass-matrix for the light neutrinos takes the following form, 

/rUee rUefj. m,er\ fA Bi 



V 



7 



V 



Ci D 
C2) 



(3.49) 
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which is generally symmetric and complex. The matrix ()3.49p can be diagonalized by a unitary 
matrix V , 



V'^ MyV = Dy = diag(m^, 771,2, 



Ml ' 



or, My = V*DyV^ , 



(3.50) 



where the diagonal matrix Dy contains three real mass-eigenvalues for the three light neutrinos. 
The mixing matrix V can be generally expressed as a product of three unitary matrices including a 
CKM-typ41j mixing matrix U plus two diagonal rephasing matrices U' and U", 



V 



u"uu' , 



with 



U' = diag(e^'^S e''f'\ e'' 



U" = diag(e*"i, e'"2, e' 



U = U23U13U12 



^ C-sC-x '^S^X ^x^ 

\SsSa + Cs 





n 








U23 = 





Ca 


-s 




^0 


Sa 


Ca 




SxC 





U12 = 




(3.51) 

(3.52a) 
(3.52b) 

(3.52c) 



where is the Dirac CP-phase. For notational convenience, we have denoted the three neutrino 
mixing angles of the PMNS matrix as. 



712, (723, t^ia; 



), 



(3.53a) 



(like in (HUI) in the above and in the rest of this paper. Consequently we further adopt the following 
notations, 



{ss, Sa, Sx) = {sm9s, sm9a, sm9x) 



{Cs, Ca, Cx) = (COSI 



COS tin , COS ( 



(3.53b) 



The diagonal rephasing matrix U' contains three Major ana phases, but after extracting an overall 
phase factor only two combinations of them persist (such as (pi — (p^ and (j)2 — ^3), which correspond 
to the two physical Majorana phases. The other three phases in U" are associated with the flavor 
eigenbasis of light neutrinos and are needed for the consistency of diagonalizing a given mass-matrix 
My. 

From the second relation of (j3.50p we can fully reconstruct all the elements of the symmetric mass- 



matrix in terms of the rephased mass-eigenvalues (rn-i, m.2, "fh^) = (m;^e~ ''^i, 17126 ''''2, m^e 
the mixing angles {Og, 9a, 9x), the Dirac phase 5d, and the rephasing phases (which do not appear 
in physical PMNS mixing matrix). 



CgCr^Tft-^ -\- S ~r ^x^ 



2„2„ 



~2iS 



(3.54a) 



^■^Here, CKM stands for the Cabibbo-Kobayashi-Maskawa mixing matrix for the quark sector, including 3 mixing 
angles and 1 Dirac CP-phase [12| . 
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{SsCa - CgSaSxC' °) + (c^Ca + SsSaS^C* °) 1712 + S^C^m^ 

{ssSa + CgCaSxe'^^ffhi + {cgSa - SsCaSxe'^'^fm2 + c\cl.m^ 





(3.54b) 




(3.54c) 


3 


, (3.54d) 


— ~ 


, (o.o4ej 




(3.54f) 



where only two Majorana phases among (/>i 2 3 (hidden in the mass-parameters 2 3) are indepen- 
dent since an overah phase factor (say e*''^3) can be taken out of U' and simply absorbed into the 
diagonal rephasing-matrix U" . For the case with a vanishing mass-eigenvalue (such as nii = in 
our present model defined in Sec. ED, only one independent phase combination [say e^^^^~'^s^] is left. 
If the light neutrino mass-matrix M,^ respects /i — r symmetry, we should have 



(3.55) 



As discussed in Sec. 12. 11 the ^—t symmetry generally allows the relation m^^ = pm^^ (with p = ± ). 
But for p = — it can always be absorbed into the phase-matrix U" by a simple shift Q3 — )• a^ + ir [cf . 
()3.54dp - ()3.54ep ]. (This shift also makes m^^ element flip a sign, but it does not affect the realization 
of /i — T symmetry and we can simply denote —mn^^ by m^^ .) Taking the three mass-eigenvalues 
(m^, 777-2, "^3) ^ independent experimental inputs without extra tuning, we can derive six mass- 
independent equations from the general conditions in Eq. ()3.55p . 



2 



C.« Cx 



I I, 



CfiSn S aCnSxG 



iSr 



I ^ 
CgCa ~r SgSaSxG 



iS 



D 



,2 2„-i2o„ 



g-j2Q!3 



(3.56a) 
(3.56b) 
(3.56c) 

(3.56d) 

(3.56e) 
(3.56f) 



The last three conditions ()3.56dp - ()3.56fp are not independent, since up to a trivial overall factor 
(l3.56dl) - (f336eD are just the squares of ()3.56ap -( |3.56bp . respectively, and (|3.56fp can be inferred from 
(|3.56ap - ()3.56bp (cf. below), where we note 0^/0°, 90° and Ox 7^ 90° due to the requirement of 
neutrino data in Table-{TJ So we only need to solve (j3.56ap -( |3.56cp . Since the data in Table-[T] already 
enforce Og ^ 0°, 90° and 6x 7^ 90°, we can deduce two equations from (|3.56ap - ()3.56bp . 



(3.57a) 
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= 0, (3.57b) 
which together give, 

3111(03 — 02) = 0, tan^a = 003(03 — 02), Sx = . (3.58) 

The 3olution Sx = also holds Eq. (j3.56cp . Noting that tan^a ^ for 6a € [0, ^] , we find that the 
// — T symmetry results in the miique consistent solution, 

{Oa, Ox) = (45°, 0°) , a2 = "3 > (3-59) 

which automatically holds ()3.56fj) . The solution ()3.59p is also the sufficient condition of realizing fx—r 
symmetry in Mi, since we can readily deduce the condition (|3.55|) by substituting Eq. (|3.59|) into the 
general reconstruction formulas ()3.54bp - (j3.54ep . 

Finally, we note that, as shown in the fJ- — T symmetric solution (|3.59p . the solar mixing angle 9s 
is independent of the /U — r symmetry and thus can take any value in the ^ — t symmetric limit. So 
we need to go beyond the /U— r symmetry for predicting Og (cf. Sec. l6.3p . though this is not the main 
concern of our present study. 



3.2. Reconstruction with Normal Hierarchy Mass-Spectrum of Light Neutrinos 

In this subsection we apply the general reconstruction formalism in Sec. l3.ll to to our model defined 
in Sec. [2] where the light neutrino mass-spectrum exhibits the normal hierarchy (NH) pattern, < 
7712 < ^3 , and has = [cf. Eq. ()2.18ap ]. As noted in Eqs. ()2.20p - ()2.2ip . the neutrino data 
(Table-[T]) requires the small mass-ratio. 



1712 



'517121 




0.17-0.19 < 1 



(3.60) 



at 2(7 level. From the relations in Eqs. (|2.2ip - (|2.24p (with p = +), we see that the small parameter 
y is related to another small quantity r in the seesaw sector and we have y = 0{r) <^ 1 . So, 
with = and y = ^ <^ 1 for the NH mass-spectrum, we can further simplify the general 
reconstruction formula (|3.54p as. 



iifi 



nine *^°2 



7n 



m3e-^2°3 



THnC 



m 



y{csCa + SsSaSxe'^'^fe-''"'^^^ + slcl 
y{csSa - SsCaSxe'^^)^e-'^'^^^ + clcl 
ysgCx {CsCa + Ss 

r„3e-^("i+"3) \-ys^Cx{CsSa - SsCaSxe'^°)e-'^'^^3 - CaCxSxe-'^^ 
yiCsC-a ~\~ SsSaSxG- ^)(CsSa SsCaSxG ^)c 



(3.61a) 
(3.61b) 
(3.61c) 
(3.61d) 
(3.61e) 
(3.61f) 
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where the Majorana phase (p^ disappears due to rrii = . We have also absorbed an overall 
redundant Majorana phase (jy^ (from U') into the redefinition of (in U"), and the only remaining 
independent Majorana phase is (j)23 , 

where j = 1,2,3. In fact, for any minimal seesaw model the low energy physical PMNS mixing 
matrix T^'MNS only contains one independent Dirac phase do and one independent Majorana phase 
(denoted as here for NH mass-spectrum). 

With these notations we may slightly rewrite the general fJ. — T symmetric solution as follows, 

{Oa, 6'x)o = 0) , = «30> (3-63) 

where the subscript "0" denotes the leading order (LO) results from the /i— r symmetric limit. 

Next, we will explicitly identify the /i — t and CP violations from the reconstruction formula 
()3.6ip . Since the existing neutrino data support the /i— r symmetry prediction ()3.63p as an excellent 
approximation [cf. Table-[T]and Eq. (jl.ip ]. we are strongly motivated to treat the small /i— r breakings, 

5a = Oa-j, 6^ = (3.64) 

as the parameters of perturbative expansion. So we can systematically analyze these breakings under 
the expansion of {6a, 6x)- For practical analysis it is enough to keep the expansion up to the next- 
to-leading order (NLO), i.e., the linear order of {6a, 6x)- Besides {6a, 6x), the other NLO parameters 
for our current reconstruction analysis will include {6m2, 6m^) and {6a-^^, 60:2, 6a^) . Since = 
in our minimal seesaw model with NH mass-spectrum, there is no NLO correction to it. Also, we 
will explicitly show in Sec. 14. II that for our expansion for the current model, the LO mass for 777-2 
vanishes (?772o = 0) as well, so we have, y = — - = , where 7712 = 77730 + "^"^2 ~ "^"^2 • 

can introduce another small ratio z = = 0{y) . In addition, as will be generally proven in 

Sec. 16. 21 the solar angle Og {= ^12) is independent of the soft /U— r breaking and thus does not receive 
any NLO correction. Furthermore, we note that the Dirac phase 6**^° is always associated by the 
small mixing parameter Sx {= 6x) and the remaining Majorana phase e*''^23 in our present minimal 
seesaw model with NH mass-spectrum is always suppressed by the small mass-ratio y {= ^) , so 
both phases only appear at NLO and thus receive no more correction at this order of expansion. 
Finally, we can summarize all independent NLO parameters in our reconstruction analysis as follows, 

( y, z, 6a, 6x, 6ai, 60.2, ^ag ) , (3.65) 

■ 1 "^2 "^"^2 "^"^2 1 <^^3 .-/-I/ N 1 r 1 ■ 1 n 1 1 

With y = — = ~ and z = = 0{y) . Lach 01 these parameters is defined as the 

mg mg mgp TTlgo 

difference between its full value and zeroth-order value. In Sec. H] we will use our soft /7 — r and CP 
breaking model (defined in Sec. l2.2p to predict these deviations, especially the correlation between 6a 
and 6x which will be probed at the upcoming neutrino experiments. 
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Under the perturbation of (|3.65p . we expand the hght neutrino mass-matrix ()3.49p up to NLO, 
My = M^°) + + (higher orders) 



/ ml? 



,(0) 



,{0) ^(0) 

I'lJ.T 

(0) 



/ 6m^J dm^/} S'm2r\ 



+ 



6m^rT J 



(3.66) 



The LO matrix M^"^ and the NLO matrix M^^^ can be exphcitly derived from our reconstruction 
formula (j3.6ip . Let us first deduce the LO mass-matrix mI)^^ from ()3.6ip . under the r symmetric 
solution (j3.63p . 



/ (0) 

V 



(0) 



m 



(0) _(o) 
(0) 

rriTT 



/Ao Bo Bo\ 

Co Do 
\ Co) 




where we have also matched to our notation of m!)'^ in (|2.4ip . So, as expected, the zeroth order r 
symmetric mass-matrix Adj)^^ only contains an overall CP-phase e~*^°2o which has no observable effect 
at this order and is needed only for consistency of the mass-diagonalization procedure (cf. Sec. l4.ip . 
From ()3.67p . it is clear that mI)^^ has the LO mass-eigenvalues (0, 0, m^Q) , while the mixing angle 
6s is not fixed at this order. 

Then, we derive elements of the NLO matrix mI)^^ from (|3.6ip . 



6mi]} 



6A 
6Ci 

= SC2 
5Bi 
6mg^ = 8B2 
(5mW = 6D 



'^'"'fiH 



1 



2 
1 
2 

1 

V2 
1 

71 



ycge 



-i{aiO+a2o) 
-«{aio+a2()) 



1 



(3.68a) 
(3.68b) 

(3.68c) 

(3.68d) 

(3.68e) 

(3.68f) 



where we have matched to our notation of 5M^^ defined in (j2.4ip . In the above formulas we have 
input the ^— r symmetric LO parameters or relations, mi = 0, (^aO) ^xo) = (f ) 0) and = c^so 
[Eq. (13:631) ]. 

According to the definitions in ()2.42p . we can uniquely decompose the elements of 6M^^ in ()3.68p 
as the IX — T symmetric and anti-symmetric parts, 5M^^ = 6M^ + 6M^ , with their elements given 
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by, 



V2 



A R 
OJDi 


_L X R 




o 
Z 




A R„ 




2 




+ SC2 




2 




-<5C2 



(3.69) 



6Ca ^ ^ = ^m3oe-'2"2o [25^ _ i(5a2 - Jag)] , 
where we may further decompose 5mI)^^ according to (|2.45p . 



(3.70) 



In the next section, we will apply the above reconstruction formulas (j3.66p - (j3.70p to match with 
()2.4ip - ()2.47p in our soft breaking model at the LO and NLO, respectively. From this, we can connect 
the seesaw parameters with the low energy neutrino observables and derive quantitative predictions 
of our soft breaking model. 

Before concluding this section, we note that so far we have presented our reconstruction formalism 
at the low energy scale. There are possible renormalization group (RG) running effects for the low 
energy neutrino parameters when connecting them to the model predictions at the seesaw scale. Such 
RG effects were extensively discussed in the literature |32j . The application to our present analysis is 
straightforward. Below the seesaw scale, heavy right-handed neutrinos can be integrated out from the 
effective theory and the seesaw mass-matrix Mjy for light neutrinos obey the one-loop RG equation 
(RGE) [32], 

dM^ 1 



dt 167r2 



aM, + c (M,(y/y£) + (y/y,)^M,jJ , (3.71) 

where t = ln(^/^Q) with jj the renormalization scale, C = — | for the SM, and denotes the 
diagonal lepton Yukawa coupling matrix in the lepton mass-eigenbasis. For the SM, the coupling 
parameter a consists of 

- -35i + 6y2 + A, (3.72) 

where {§2, Ut, A) stand for the SU {2)l weak gauge coupling, the top Yukawa coupling and Higgs self- 
coupling, respectively, which are all functions of t at loop-level. Thus, the RGE for mass-eigenvalues 
rrij {j = 1, 2, 3) can be derived as [32] . 

dm, 1 



[a + 2CylTj]mj , (3.73) 



dt 167r2 

= 2s1sl- 'iSxSsCsSaCa cos 5d + 2sIcIcI, 
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J"2 = 2c^S^ + ASxSsCsSaCa COS 6d + 2s^S^C^ , 
^3 = 2c2c2, 

where = O(10~^) is the largest lepton Yukawa couphng. We see that the flavor-dependent term 
on the RHS of (|3.73p is suppressed by = O(10~^) relative to the universal a-term, and thus 
completely negligible for the present study. Hence, from the RGE ()3.73p we can express the running 
of the mass-parameter rrij from the scale to n , to good accuracy, 

iT^iM = x(m,/^o)'^7(^o) > (3.74a) 



^o) - exp 



1 



167r2 



a{t') dt' 







(3.74b) 



where t = ln(/i//iQ) . For current study we will set, (/1q, /i) = {Mz, Mi) , where Z boson mass Mz 
lies at the weak scale and the heavy mass Mi represents the seesaw scale. 

For the minimal neutrino seesaw with NH mass-spectrum, = < m2 <^ , we note that 
the zero-eigenvalue and the ratio y = ^ do not depend on the RG running scale /.i . So we 
can derive the running of the two nonzero mass-parameters from weak scale to seesaw scale, 

m2 = m2{Mi) = Xi'm2{Mz) = yfh^, (3.75a) 
mg = m^{Mi) = Xi'^siMz), (3.75b) 
Xi = x(Mi,Mz), (3.75c) 

and accordingly, A^(Mi) = Xi^s{Mz) and Aa(Mi) = Xi^a{Mz) ■ The RG running factor Xi = 
x{Mi, Mz) can be evaluated numerically, and depends on the inputs of initial values for the weak 
gauge coupling a2 = (7|/(47r) , the top-quark Yukawa coupling and the Higgs boson mass Mh , 
via the combination a in Eq. (j3.72p . The existing electroweak data gives, ^(M^) = 29.57 it 0.02 , 
mt = 173.1 ± 1.4 GeV, and the Higgs mass range 115 ^ ^ 149 GeV for the SM at 90%C.L. 
[13\ I33j . So the universal running factor xiMi, Mz) is found to be around 1.3 — 1.4 for Mi = 
10^^ - 10^^ GeV and Mh = 115 - 149 GeV. In Sec. HE] this running effect will be numerically 
computed for our analyses. Other running effects due to the leptonic mixing angles and CP-phases 
are all negligible for the present study since their RGEs contain only flavor-dependent terms and are 
all suppressed by = 0(10""^) at least [32] . 

For the analyses in Sec.|l][5l we will first run the low energy neutrino parameters in the cur- 
rent reconstruction formalism up to the seesaw scale Mi , and then match them with our model- 
predictions. The relevant places for including such RG effects are just to replace the light neutrino 
mass-eigenvalues (?7i2, m^) at the low scale by (rrig, fh.^) at the seesaw scale Mi . 
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4. Predictions of Neutrino Seesaw with Common Soft fJ^ — r 
and CP Breaking 



In this section we systematically derive the predictions of our common soft /i — r and CP breaking 
model (Sec l2.2p for the low energy neutrino observables by using the reconstruction formalism (in- 
cluding the RG running effects) in Sec. El Especially, we will analyze the correlation between the two 
small fi — T breaking parameters 6x (= ^13 — 0) and 6a (= O23 — j) ■ Other correlations of 6x with 
the Jarlskog invariant J or the Qv(3(3 decays parameter \m^^\ are also analyzed. 

4.1. Correlation between 6^^ and 023 — 45° 

In this subsection, we first analyze a general feature of the NLO mass-matrix which predicts 

5x oc 5a and imposes a lower hound on 5x for any nonzero 5a , independent of all details of the /x— r 
and CP breakings. 

Inspecting the reconstruction formula ()3.69p . we see that 6x and 5a are only contained in the 
/X — r anti-symmetric parts 6Ba and 5Ca, respectively. Using these two formulas we compute their 
ratio, 

5Ba V24e-^("io-"2o+'5i5) 

— - = - — . (4.76) 

5Ca 25a - i{5a^ - 5a^) 

On the other hand, from Eq. ()2.48p [cf. also Eq. ()2.35p ]. our soft breaking model predicts this ratio 
to be real, 

5Ba a 1 , , „„, 

5Ca b + c k 

which does not depend on any detail of the /i — t and CP breakings. According to our general proof 
(cf. Sec. 16.21 for detail), this ratio k is related to the solar mixing angle Og via Eq. ()6.15ip . So we can 
rewrite (0777]), 

g = ^'a-".. (4.78) 
where = it denotes the sign of k . So, from (j4.77p and (|4.78p we deduce two conditions, 

5a = 5x coi 6s cos{aiQ- a2Q + 5 D)Pki (4.79a) 

5a2 — 5a.^ = 2tan(a;^o ~ ^^20 + fe) (^a ■ (4.79b) 

Eq. (j4.79ap shows that at this order the two small /i— r breaking parameters are proportional to each 
other, 5x c<. 5a ■ Strikingly, due to | cos(ci]^g 0^20 ~^ ^d)\ ^ ' can further derive from (j4.79ap a 
generic lower hound on 5x , for any nonzero 5a, 

5x ^ 1 5a I tan 0,, (4.80) 

where 5x = ^13 G [0, f ] in our convention. 
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We stress that the above derivation of ()4.79ap -( )4.80p depends only on the ratio of two ^— r anti- 
symmetric elements 6Ba/5Ca in (j4.76|) . and our model predicts this ratio as 1/k in (|4.77|) . which 
originates from the Dirac mass- matrix rrij-, alone and does not depend on any detail of the r and 
CP breakings. Hence, Eqs. (|4.79ap and (|4.80p reflect the general feature of our model. 

4.2. Full NLO Analysis for Low Energy Neutrino Observables 

In our model, both /z — r and CP violations arise from a common origin in the seesaw Lagrangian, 
characterized by the unique soft breaking parameter Qe^^ (appearing at the NLO of our perturbative 
analysis). Thus, the small /i — r breaking parameters {5a, 5x) as well as the low energy Dirac and 
Majorana CP phases {6d, 4>23) the light neutrino mass-matrix are controlled by ^ ^-nd w . In this 
subsection, using the general reconstruction formalism (Sec. [3]) for diagonalizing the light neutrino 
mass-matrix at the NLO, we will systematically derive predictions of our model for these low energy 
observables and their correlations. 

Let us first have an examination on the reconstructed LO mass-matrix mIP^ in (j3.67p . Comparing 
the ()3.67p with our model prediction ()2.40ap at the same order, we can deduce, 

where we have evolved the relevant low energy mass-parameter in (j3.67p up to its value at the seesaw 
scale via ()3.75p . the symbol = ± denotes the sign of r as before, and Q3Q = due to the fi — r 
symmetric solution (|3.59p . The fact that m.^Q is the largest mass-eigenvalue at the LO shows that 
light neutrinos pose NH mass-spectrum in our model. Also, in this expansion the solar mixing angle 
6s is not fixed at the zeroth order, and it will be determined at the NLO below. 

We then inspect the light neutrino mass-matrix 5M^^ at the NLO, as given by our soft breaking 
model in Eqs. ()2.42p - ()2.47p and by our reconstruction formulas in Eqs. (|3.68p - (|3.70p . We can match 
the two sets of formulas for the /x — r symmetric elements, 

„2 

5A = -^r = m3oys?e-^2("io+'^23), (4.82a) 



Ml, 







SBs = = --1-^30 y.,c.e-('^io+^.o+2*.3), (4.82b) 

5Es = = la3oyc^-2('^-+'^-), (4.82c) 

^E's = .^tfS. r = ^m3oe-2^- [z - z{6a, + 6a,)] , (4.82d) 



4(2-X)2Mio 
and for /x— r anti-symmetric elements 



= -ii^ - = . (4.83a) 
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SCa = - ^l^-XW.o ^ ^"^3oe-*'^^«[2'^a-i(fc2-'^«3)], (4.83b) 

where we have run the low energy mass-parameter up to fh^Q at the seesaw scale M\ via 

Eg. (1X75]) . 

We first analyze the ^ — t symmetric Eqs. ()4.82ap -( |4.82d|) . Inspecting the model predictions for 
{5A, 6Bs, 6Es) on the left-hand-sides (LHS's) of (l4.82al) - (H:B2cD . we observe that {6A, 6Bs, 5Es) 
are all real. This requires the RHS's of (|4.82ap -( )4.82cp to have vanishing imaginary parts. We thus 
derive, 

2(aio + 4>23) = f^iT^ ' "10 + "20 + 2</'23 = ^^2^ > (4-84) 
where u-^tt, ngvr = 0,7r (mod 27r). From these we have, 

"io + <?^23 = -y-' "10 -"20 = (4-85) 

where nvr = (nj^ — 712)71 = 0, vr (mod 27r). Then, we compute the ratio of ()4.82ap and (]4.82bp . and 
arrive at 

tan^, = (-)-+i^, (4.86) 
b + c 

where we can always ensure tan^^ > 0, i.e., we set n = 1 for a/{h + c) > 0, and n = for 
a/ (b+c) < . We note that due to the NH mass-spectrum, we have a small parameter r = ©(rrig/w-g) 
in our present model [cf. ()2.18ap with p = +] and a corresponding small parameter y = "012/ in 
our reconstruction formalism [cf. ()3.60p ]. As a result of expanding r and y in the NH scheme, it is 
clear that the solar mixing angle 9s is not fixed from the LO mass- matrix M^f^ in (|2.4Uap or (|3.67p . 
but determined at the NLO as above. (In Sec. l6.H[6.2l we will generally prove that the Og formula 
(j4.86p actually does not depend on the inclusion of /x — r and CP violations at the NLO.) 

Comparing ()4.86p with the tan^^ in ()4.77p - ()4.78p . we can fix the sign of A; as = {—)^~^^ ■ 
With the second relation of ()4.85p . we can further simplify ()4.79ap - ()4.79bp as 

-r- = — cot 9s cos 6d , (4.87a) 

Ox 



-'X 



— Sa^, = 2 tan 5z) (5a, (4.87b) 

where the sign of ()4.87bp does not depend on the choice of n . Since our convention already ensures 
cot 0s and dx (= sin^^:) to be both positive, the sign of the deviation 5a just equals the sign of 
— cos 6d ■ 

With Eqs. KSTEh and (jTOD - (jT85D . we can further solve from Eqs. (Ii:82ip - (|4.82dp . 

= (-)"is^m3oMioF, (4.88a) 
a{b + c) = (-)"2+^\/2s,c,m3oMior, (4.88b) 
{b + cf = (-)"i2c2m3oMioy, (4.88c) 
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i|(2-X)Mio|m30, 



(4. 



where for convenience we have defined a ratio, 



Y 



y 



(4.89) 



bmce MiqY = 2/M22 > and a^, (6 + c)^ > , we deduce, n^^ = , which further leads to n2 
[cf. the comment below ()4.85p ]. With these we can solve the expressions for (a, 6, c) , 



-n 



h + c = pj^Csyj2m^QMioY , 
6-c = p_wJ|(2-X)Mio|m3o , 



where p-t = it denotes the signs of b± c. Combining Eqs. ()4.8ip and (|4.82d[) . we derive, 

———r = z - i {6a2 + da^) , 
which can be decomposed into real and imaginary parts. 



\X\ 



6ar, + 6ao 



4|2-X|2 



412-X 



(2cos2w — \X\ coscj) r , 
Y (2 sin 2(jj — \X\ sin a;) r . 



(4.90a) 
(4.90b) 

(4.90c) 

(4.91) 

(4.92a) 
(4.92b) 



Next, we analyze the fi—r anti-symmetric equations ()4.83ap -( )4.83bp for 6mI)^\ With Eqs. ()4.81bp . 
K90\i and ()OTD - (|T85D . we can rewrite (|i:83all - (l4.83bD as 



X^Y 



2-X 



-P+P- 



X^Y 



2-X 



25a — i{Sa2 — 5a^) . 



Taking the ratio of the two sides of (j4.93ap - ()4.93bp and making use of ()4.87bp . we deduce 



= (-)"-"2+i cote, cos 

Ox 



(4.93a) 
(4.93b) 

(4.94) 



which is consistent with Eq. (|4.87ap due to n2 = —n as we derived below Eq. ()4.88p . Then, we may 
rewrite ()4.85p as 



a 



10 



a 



20 



nvr 



^23 ' 



(4.95) 



where we have used rii = derived below (j4.88p . With this we can derive, from ()4.81bp . the 
Majorana phase angle (^'23 in terms of the original soft CP-breaking phase angle u , 

C sincj 



tan 2023 



2r — cos uj 



(4.96) 



27 



Then, from Eq. ()4.93p we finally solve, 



Sa 



27r 

1 

2 



2-X 



s.r 



2-X 
X^Y ^ 



cos 6d Cs\r\ , 



2-X 



smSo Cs\r\ , 



(4.97a) 
(4.97b) 
(4.97c) 
(4.97d) 



where we have defined the phase angles w, 6d £ [0, 27r) and the mixing angles {6s, Oa, Ox) G [0, f ], 
which requires the sign product p^p^p^ = + for consistency. From (j4.97p we can reproduce the 
relations (I4.87ap - (l4.87bl) in Sec. 14.11 Furthermore, using ()4.97bp - ()4.97cp we can explicitly express 5x 
and 5a in terms of r and C, , 

g ^ y/ySsC 

2 [C^ - 4rC cos 5d + 4r2] 4 

. _ -^CsCos6dC 



(4. 

(4.98b) 



2 [C^ - 4rC cos 6d + 4r2] 4 

where y = ^ ~ y^A^/A^, G [0, 2tt) and < C < 1 in our convention. Since — 4rCcos5i:) + 
4r^ = (2r — Ccos5/))^ + {Csindo)'^ > , so the denominators in (j4.98p are well-defined. 

It is useful to note that reversing (j4.98ap we may resolve the seesaw parameter r as a function of 
C and 5 d , 



cos 5d 



16 St 



sin^ 6d 



(4.99) 



where the parameters Ss (= sin 6*12) , V [ = ^ — \J ^s/^a j and 5^ (= ^13) have been measured 
by the oscillation experiments (Table-[T|). The square-root in ()4.99p has to be real, so we have the 
condition. 



452 

sm^ul , 



(4.100) 



which will put a lower bound on for nonzero bx and sin 5/) . 

Inputting the neutrino data for Qg and (A^, A^) in Table-[T] and scanning the phase angle bo G 
[0°, 360°) , we can plot (= 6x) and 6*23 - 45° {= 6a), from (|4.98ap - (|4.98bp . as functions of the 
parameters C and 6r) . This is shown in Fig. [1] with the experimental inputs varied within 90% C.L. 
ranges and with the natural regions of \r\,(^ € [0, 0.6] . Inspecting Fig.[T^ and Fig.[T]:, we see that for 
C G [0, 0.6] , an upper bound ^13 < 6° holds for most of the parameter space. Fig.[T]3 depicts how 
the deviation 623 — 45° behaves as a function of C via Eq. ()4.98bp in which all measured observables 
are constrained by the 90% C.L. data and the variable ( needs to ensure that 613 {6x) obeys the 
current oscillation limits (Table-[T]) via Eq. ()4.98ap . As a result, we see that both Fig.[T^ and Fig. [1)3 
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Figure 1: The predictions of ^13 and ^23 — 45° as functions of the soft jJL—T breaking parameter C, and 
CP-phase angle Sd, are shown in the upper plots (a)-(b) and lower plots (c)-(d), respectively. The 
experimental inputs are scanned within 90% C.L. ranges and the Dirac phase angle 6d G [0°, 360°) , 
with 2000 samples. The shaded region (yellow) denotes the current bounds (90% C.L.) on ^23 — 45° 
according to Table-dl 



confined our parameter space into C ^ 0-1 region. We further note that in Fig.[T]i there is a positive 
peak-region of ^23 — 45° for the parameter space 120° ^ (5_d < 220° ; but Fig.[T]: requires that around 
the same range of 5d the ^13 is much suppressed, and the two peak-regions of ^13 in Fig.[lJ; around 
(5d ~ 50° and ~ 310° just correspond to the negative deviation, ^23 — 45° < . 

With (j4.98ap -( )4.98bp we can plot the correlation between the two ^—t breaking mixing angles ^13 
and ^23 — 45° . This is shown in Fig. [21 where we have varied the measured parameters within their 
90% C.L. ranges and the Dirac-phase angle 5d G [0°, 360°) , as well as C;'^ ^ [0) 0-6] . The current 
90% C.L. limits on ^13 are shown by the shaded region (yellow), while the ^13 sensitivities of the 
upcoming Double Chooz [T3j, RENO [16] and Daya Bay [15] experiments to are shown by the three 
horizontal (red) lines at 90% C.L., as 5.0°, 4.1° and 2.9° (from top to bottom), based on three years 
of data-taking. From Fig. [21 we see that our model generally predicts ^13 < 6° ; the Double Chooz 
experiment [T3] could marginally probe the predicted ^13 and the RENO experiment [16] can be a 
better job. The NOi^A experiment [T7] is also expected to probe ^13, with a sensitivity of ^13 ~ 4.4° , 
slightly lower than RENO. It is clear that the Daya Bay experiment 115^ holds the best sensitivity 
and can probe the large central region of our predicted parameter space. (For more details of the 
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Figure 2: The correlation between ^13 and ^23 — 45°, based on Eqs. ()4.98ap -( )4.98bp . wliere we have 
scanned the experimental inputs within 90% C.L. ranges and the Dirac-phase angle 6d S [0°, 360°) , 
with 2000 samples. The sensitivities of Double Chooz [2], RENO jl6j and Daya Bay jl5j experiments 
to are shown by the three solid horizontal lines at 90%C.L., as 5.0°, 4.1° and 2.9° (from top 
to bottom), based on three years of data-taking; the dashed horizontal line gives Daya Bay's future 
sensitivity (2.15°) after six years of running j42j. 

experimental sensitivities to ^13 , see a nice comparative analysis in |42j.) 

We note that the sharp edges on the two sides of the allowed parameter space in Fig. [2] are 
essentially determined by the lower bound given in ()4.80p . 6x ^ |5a|tan0s; where the measured 
parameter tan^s ~ 0.67 (Table-{T|) just corresponds to the slopes of the sharp edges which are nearly 
straight lines. This means that for any measured value of 023 — 45° 7^ , the Fig.\^ imposes a lower 
hound on 612,, which will be probed by the ^13 experiments such as Daya Bay and RENO. This is 
a really interesting and encouraging prediction for the upcoming neutrino oscillation experiments 
which will probe the /i — r violating observables ^13 — 0° and ^23 — 45° to much higher precision. 
Note that the current oscillation data already favor the central values of ^23 to be smaller than 45° 
(Table-[T|) and this feature is quite robust (cf. footnote-l2] and Ref. [7]). Hence, our findings strongly 
encourage the experimental efforts (such as MINOS 120^ and T2K 118^) to further probe 623 as precise 
as possible, hopefully to the similar level of accuracy as the present 9i2 measurement. 

Next, we analyze our model predictions for the low energy CP-violation (via Jarlskog invariant) 
and the neutrinoless double beta decays (via the element \mee\ of M^). In our construction (Sec. l2.2( ). 
the original soft breaking CP-phase e*" is the source of both low energy Dirac and Majorana CP- 
violations via the phase angles 6 d and 4*2^ ■ 
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Figure 3: Correlations of ^13 (in degree) with the Jarlskog invariant J [plot-(a)] and with the 
neutrinoless-doublet-beta-decay observable Mg^ [plot-(b)]. Each plot has computed 1500 samples. 

The Dirac CP-phase 6^ will manifest itself in Jarlskog invariant J while both 6d and the Majorana 
CP-phase (j)2^ will appear in the neutrinoless double beta decay observable Im^J . The Jarlskog 
invariant J can be written as 



J = — sin 20<i sin sin 20a; cos 20a; sin 



sin 20, sm5D + 0{5l,6l), 



(4.101) 



where as defined earlier, 6x = 0x and 5a = 9a — j ■ 

Then we analyze the neutrinoless double-beta decay. Our model predicts the NH mass-spectrum 
with = 0, so from (|3.61ap we can derive the mass-matrix element |mgg| for neutrinoless double 
beta decays, 



mo 



IE*'. 



2 



y •^s'-x^ 



- rn,Jy^si + 2y5lsj cos 2(60-^3) + i^x " '^^jy^^ 



(4.102) 



From Eqs. (|4.95p . ()4.97ap and noting ago ~ '^30 ' derive the phase-angle combination, 

i^D - (t>23) = ("30 - w) + (2 + n)7r . 
Thus, we can compute, by using Eq. (|4.81bp . 

cos2(5d-023) 



(4.103) 



-iluj 



{2-X) 



\2-X\ 
2r cos 2(5/5 — C cos 



(4.104) 



y/C^-ArC cos 6d +4r2 
where we have used the solution lo = 2-k — in (j4.97ap . 

In Fig. [3^ we show the correlation between 0i3 and the Jarlskog invariant J, and in Fig. [3)3 the 
neutrinoless doublet beta decay observable Mge (= jw-ggl) depicted. In plotting this figure, we have 
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Figure 4: The mass ratio M2/M1 of right-handed heavy neutrinos as a function of the soft /i — r 
breaking parameter C, , with the experimental inputs varied within 90% C.L. ranges [plot-(a)] and 2(j 
ranges [plot-(b)], respectively. Each plot has computed 2000 samples. 



used Eq. ()4.98ap where ^13 is taken as a function of r and C,. All other measured parameters in the 
formulas of ^13, J and Mee are scanned within their 90%C.L. ranges, while the Dirac CP-phase 5d is 
varied in [0, 2tt). The shaded region (yellow) in this figure is allowed by the current data at 90% C.L. 
From Fig.[3)[a), we see that a nonzero J can place a lower bound on ^13 in our model. This is because 
I sin(5£)| ^ 1 and thus we can deduce from Eq. (|4.10ip . 

4|J| 



sin 20O 



(4.105) 



Fig-E^b) shows an upper bound Mge ^ 3.4 meV which is quite small (as commonly expected for 
schemes with NH mass-spectrum) and poses a challenge to the future neutrinoless double beta decay 
experiments j31j . 



4.3. Analysis of Seesaw Parameter Space 

In this subsection, we further analyze the allowed seesaw parameter space, which will also be needed 
for our leptogenesis analysis in Sec.[5l 

The mass-eigenvalues of the right-handed neutrinos can be directly derived from diagonalizing 
the Hermitian matrix M^^Mj^ or Mj^M^^, under the expansion of ( and r ^ 



Ml 
M2 



M22 
2M22 



— rCcosui + -C'^ 
4 



1 — -(2r -|- <^cosw) 



So we derive the mass ratio, 



M2 



A - {2r + ( cos uj) 



Ml 



1 ^ (2r — ( cos Ci;)^ -|- sin^ 



(4.106a) 
(4.106b) 

(4.107) 



^*We can also diagonalize the complex matrix M_r by using the unitary transformation matrix Vr , as shown in 
Appendix- A. 
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As shown in Eq. ()4.99p . the seesaw parameter r can be resolved as a function of ^ and uj with 
(ss, y, 6x) from the oscihation data. In Fig. U] we plot the ratio M2/M1 as a function of the ^ — t 
breaking parameter C, , where we have scanned the parameter space a; G [0, 2tt) . As we noted earlier 
in Eq. ()2.23ap . the neutrino data already favors M2/M1 ^ 1 . For our leptogenesis study in Sec. [5] 
we will also require M2/M1 ^ 5 (as indicated by the (red) horizontal line in Fig.[l|), so that lepton 
asymmetry is mainly produced by the decays of the lighter right-handed Majorana neutrino with 
mass Ml. 
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Figure 5: Seesaw scale Mi as a function of the elements |fe ± c| in the Dirac mass-matrix m^, is 
shown in plot- (a) and (b), where the shaded regions correspond to the natural perturbative region 
|6 ± c| G [1, 300] GeV. This puts an upper bound, Mi < 1.4 x 10^^ GeV, from plot-(b). The natural 
parameter-space in the Mi — M2 plane is depicted in plot-(c), where the shaded region is bounded 
from the right by the upper limit of Mi as inferred from plot-(b) and from the left by the lower limit 
of Ml as imposed by (j5.127p via the leptogenesis in Sec. 15.21 



From Eqs. (gSHc]), (|4.88dp . (|4.1()e)a|) and (I4.97ap . we can connect the seesaw scale Mi to the 
element of Dirac mass- matrix mp, , 

1 

(4.108a) 



(4.108b) 

where the Dirac mass-parameters, 6ibc = {y^^y^v / , arise from the Yukawa interactions (with 



Ml 



+ c 



— rC, cos 5d + 
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couplings and for b and c, respectively), and the seesaw parameter r is given by Eq. ()4.99p . 
Thus, we can plot Mi as a function of the magnitude of the Dirac mass-parameters |6 it c| in 
Figs. [5^ and[5]3, where we have input the measured quantities in their 90%C.L. ranges and scanned 
the allowed parameter space for (^, 5d), with 1200 samples. 

We note that the Yukawa couplings y^, and cannot be too small (to avoid excessive fine-tuning) 
or too large (to keep valid perturbation). So, we will take the Dirac mass-parameters |6 ± c| in the 
natural range [1, 300] GeV, corresponding to the Yukawa coupling y^ 's no smaller than 0(10"^) and 
no larger than 0{y^), where y^ = ^plm^^jv ~ 1 is the top-quark Yukawa coupling in the SM. This 
natural perturbative range of |6 ± c| is indicated by the shaded area in Figs.[5^-b, which results in 
an upper limit on the seesaw scale M\ due to the perturbativity requirement. We see that Fig. [5)3 
puts a much stronger bound on M\ than Fig. [5^ because Eq. ()4.108ap has much larger uncertainties 
due to scanning the parameter-space of C and . We may also resolve M\ from Eq. (|4.88ap . so 
M\ is expressed in terms of the Dirac mass-parameter a , 



Ml 



2 - 



— r( cos 6d + 



(4.109) 



'2 

which has similar uncertainties to Eq. (|4.108ap and thus does not provide better constraint on Mi . 

Using the eigenvalue formulas (j4.106ap - ()4.106bp . we have further plotted the viable parameter- 
space in the Mi — M2 plane, as shown in Fig.[5h with 1200 samples, where the shaded region is 
bounded from the right by the upper limit of Mi as inferred from Fig. [5)3 and from the left by the 
lower limit of Mi as imposed by ()5.127p via leptogenesis (which will be derived later in Sec. l5.2p . 
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5. Origin of Matter from Soft fi — r and CP Breaking 



In this section, we derive the predictions of our soft breaking seesaw model for cosmological matter- 
antimatter asymmetry (baryon asymmetry) via thermal leptogenesis |24^ [25] . We fully reconstruct the 
leptogenesis CP-asymmetry from the low energy Dirac CP phase and establish the direct link between 
the cosmological CP-violation and the low energy Jarlskog invariant. We analyze the correlations 
of the leptogenesis scale with the low energy observables such as the leptonic Jarlskog-invariant J 
[27] and neutrinoless double-beta decay parameter Mge [31]. We also deduce a lower bound on the 
leptogenesis scale for producing the observed matter-antimatter asymmetry. 



5.1. Matter-Antimatter Asymmetry via Leptogenesis in Neutrino Seesaw 

The universe is populated exclusively with matter instead of antimatter, and the baryon density 
( ^ ~ ) relative to the photon density is found to be small but nonzero [26] , 

r/o = — ^ = (6.21 ± 0.16) X 10~^°. (5.110) 

Generating a net baryon asymmetry requires three conditions a la Sakharov p3] : (i) baryon number 
violating interactions, (ii) C and CP violations, and (iii) departure from the thermal equilibrium. 
The standard model (SM) could not provide the observed baryon asymmetry due to having too 
small CP-violations from the CKM matrix and the lack of sufficiently strong first-order electroweak 
phase transition [35j . But, in the seesaw extension of the SM, the thermal leptogenesis [23] has CP- 
violation phases arisen from the neutrino sector and the lepton number asymmetry produced during 
out-of-equilibrium decays of heavy Majorana neutrino Nj into the lepton-Higgs pair £H and its CP- 
conjugate iH*. Since the nonperturbative electroweak sphaleron [36] interactions violate B + L |37| 
but preserve B — L, the lepton asymmetry is partially converted to a baryon asymmetry [38l |39], 

Vb = jNLl = -y<> (5-111) 

where ^ is the fraction of B — L asymmetry converted to baryon asymmetry via sphaleron pro- 
cess [38], C = i8NF + 4:NH)/{22NF + l3NH) with Np {Nh) being the number of fermion gener- 
ations (Higgs doublets). The SM has {Np, Nh) = (3, 1) and thus ^ = 28/79. The parameter 
/ = N^'^^/N* = 2387/86 is the dilution factor calculated by assuming standard photon production 
from the onset of leptogenesis till recombination |39j . The contribution from decays of the heavier 
right-handed neutrino (A'^2) will be washed out through the thermal equilibrium, only the lightest 
one (A'^i) contributes effectively to the net lepton asymmetry so long as Mi <^ M2. This is prac- 
tically realized by requiring M2/M1 ^ 5 as indicated by the (red) horizontal line in Fig.Ul Note 
that M2/M1 ^ 1 is also consistent with the condition ()2.23ap we derived in Sec. 12.11 Then, the net 
contribution to lepton asymmetry n[ is given by [39] , 

n[ = ^KfCi, (5.112) 
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and thus we have the final baryon asymmetry, 

3e 



Vb 



4/ 



Kre 



(5.113) 



with d = 3^/(4/) ~ 0.96 X 10^2 ^ 

Kj- is an efficiency factor that depends on how much out-of- 
equilibrium A'^i-decays are. It is deduced from numerically solving the Boltzmann equation |39t I4U| . 
Practically, one can derive useful analytical formulas for Kj by fitting the numerical solution of the 
Boltzmann equation. For convenience, we will use the following fitting formula for [50] ^ 

^■^^ .-^ .-^ y 1 n-3 pV 

(5.114) 



m 



0.55 X 10^3 eV 



+ 



TTli 



where 



mi 



(mj^m 



Ml 



(5.115) 



and rrij^ = m^Vn , with Vr determined from the diagonalization of Af^ (cf. Appendix-A). The CP 
asymmetry parameter can be expressed as 



r[iVi ^ £H] -T[Ni^ m* 



D"''D)12\ 



(5.116) 



T[Ni^m] + T[Ni^m*\ ^T^v' \MiJ (ml^m^yn 

where v is the vacuum expectation value of the SM Higgs boson. Note that in the mass-eigenbasis of 
(A'^i, A'^2) their Yukawa couplings and the corresponding effective Dirac mass-matrix fhj^ (= m^Vn) 
violate CP due to the complex Mr and thus its diagonalization matrix Vr in our soft breaking model, 
even though the original weak-basis Yukawa couplings and nij^ are real. It is this complex that 
differs the decay width T[Ni IH] from T[Ni IH*] . For the SM, the function F{x) in (j5.116p 
takes the form, 



Fix) 



X l-(l + x^)ln 



1+x^ 



+ 



for X ^ 1 . 



(5.117a) 



(5.117b) 



As explained earlier, in our numerical analysis we will require the mass ratio M2/M1 ^5, so we 
see that ()5.117bp is a rather accurate approximation to (j5.117ap . 



5.2. Model Predictions for Leptogenesis 



In the formulas (|5.114p and (|5.116p . we have defined a mass matrix which is connected to the 
original Dirac mass matrix , 



rriD 



rriDVR 



e*Ti [c^b + Sj^e'^c] e*^2 [cj^e^^c ■ 



(5.118) 



Other fitting formulas to tlie exact solution of also exist in the literature [39] and the resulted values are 
quite close in the relevant range of mi . 
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where the unitary matrix Vr is given by (|A.170a|) in Appendix- A from the diagonaUzation of M/j. 
With this we compute, 



D)l2 



(a^ + 26c)(4e^^-4e-'^)e 



-j(tJ+7l-72) 



(5.119a) 



(5.119b) 



Then, we can derive the parameter nii in (|5.115p . 

(6 - c)2 



mi 



as weh as the imaginary part. 



Ml 



^2 



\Y\ - -\2-X\ 



sincj (4r coscj — C) 
\/ C"^ — 4rC cos uj + 



(5.120) 



(5.121) 



With the data of Table-[T] and using Eq. ()5.120p , we see that the light neutrino mass-parameter nii 
falls into the range, 0.045 < m-^ < 0.053 eV, at 3a level. This means that in the formula ()5.114p the 
first term dominates and the second term is negligible. 

Then, we further derive the CP-asymmetry parameter e^^ from (|5.116p . 



m- 



5M1 3(4y - v^C^-4rCcos(5D-F4r2 



(4r cos 6d — C) sin ^dC^ , 



(5.122) 



^ Attv"^ 128(C2-4rCcos5z)-F4r2) 

where the solution (|4.97ap . oj = 27r — So , is used to replace uj hy 6d ■ As expected, Eq. (|5.122p 
proves that in our model the leptogenesis CP-asymmetry is completely reconstructed from the the 
low energy Dirac CP phase 5£, . 

Combining the above formulas (|5.113p . ()5.116p . ()5.117p and ()5.122p . we finally derive, 

2 

(4r cos 5z)-C) sin (5DC^ (5.123) 



J}b_ 
Ml 



ITLn 



4y - \/C^ - 4rC cos -h4r 2 



^^4^ 128(C2-4rCcos SD+^r"^) 



where we note (^"^—ArC, cos dD+^r"^ = {2r—( cos 5d)^ + sin^do > . As WMAP data (|5.110p requires 
the baryon asymmetry r]^ to be positive, we have the constraint, 



(4r cos 5/) — C) sin (5/5 > 0. 



(5.124) 



Substituting the expression (|4.99p into the positivity conditions (|5.124p . we derive the following the 
allowed regions of 6d with possible constraints after a systematical analysis (cf. Appendix-B), 



for r 



.„.,0,^)u(-,)u(|.f); 



(5.125a) 



for r = r_|_ 7^ r_ , 
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Figure 6: Pictorial summary of all constraints (j5.125ap - (j5.125cp from the positivity condition ()5.124p 
is depicted in each octant of the CP-phase bn G [0, 27r) , where the octants checked by / au- 
tomatically hold this condition, the octants checked by x have no solution, the octants checked 



by > subjects to the condition cos (5d > 
cos^(5d < 44^ 



and those checked by < are constrained by 



s4 y2(^2 ■ The radial solid-lines in each circle denote the corresponding 5d values which 
automatically hold (j5.124p and the radial dashed- lines correspond to the 5d values excluded by the 
positivity condition (j5.124p . 
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5n e (0, -J 
for r = r_ ^ r+ 



^•l' 2)>^ T'^ 



Svr 

Ulvr,- 



U 



T' T 



Svr 7vr\ 



vr 



vr Svr 
2' T 



U 



Svr \ / Svr Svr \ 



"Svr 


7vr' 














T' 


T 





(5.125b) 



(5.125c) 



where r_|_ (r_) corresponds to taking the plus (minus) sign in the formula ()4.99p . Also, the regions 

4 St 



with subscript ">" require imposing the condition cos^ 5d > 4~fe^ ) and those with subscript "<" 



are constrained by the condition cos^ < 4[ ttpi ■ 

The positivity solution (|5.125p is pictorially summarized in Fig.O We see that the values 5d = f 
and 5/5 = vr are disallowed in all three cases (represented by dashed lines in Fig. [6]). This feature 
is reflected in Fig. [7^, in which a small region around (J^ ~ ^ or (5/) ~ vr is excluded. The similar 
feature will be exhibited more clearly by Fig. [11] in Sec. 15. 21 later. 

Using Eq. (|4.99p we can eliminate r in terms of (C, w) • Finally we can plot, in Fig. [71 the ratio 
i]g/Mi as a function of Dirac CP-phase 5d via Eq. (j5.12Sp . where all experimentally measured 
quantities are scanned within their 90% C.L. range, for 2000 samples. We see that Fig. [7] shows a 
robust upper bound, 



Ml 



< 2 X 10"^^ GeV 



(5.126) 



Given the observed value of rj^ in (jS.llOp . we infer a lower bound for the leptogenesis scale Mi , 

Ml > (2.9-3.3) X 10^3 GeV, (5.127) 
where the experimental value of rj^ is varied within its 90% C.L. range. 
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Figure 7: Ratio rj^/Mi is shown as a function of soft breaking parameter C in plot- (a) and as a 
function of Dirac CP-phase 6d in plot-(b), respectively, where all experimental inputs are scanned 
within their 90%C.L. ranges, for 2000 samples. 



5.3. Direct Link of Leptogenesis with Low Energy Observables 

Since the successful leptogenesis puts an additional nontrivial constraint (|5.125p on the parameter 
space, we inspect the correlation of ^13 versus 623 — 45° in Fig. [21 We find that it is now altered as 
in Fig. [SI after we impose the condition (j5.125p as well as M2/M1 ^ 5 . It is interesting to note that 
the successful leptogenesis results in a general lower bound on the mixing angle ^13 , requiring 



^13 



> 1° , (5.128) 



even for the region around 623 = 45° . This bound is still relatively weak due to the cubic power- 
dependence of r]^ (X (^^ (X 9f^ in ()5.123p for the present NH mass-spectrum of our model, which 
means r]^ is not sensitive enough to the small mixing angle ^13 . Other extensionq^ of our minimal 
model may lower the power-dependence of rj^ on ^13 and thus enhance this lower bound. 

Then, in Fig. [9] we replot the the correlations of ^13 with the Jarlskog invariant J and the 
neutrinoless double beta decay observable M^e , respectively. This should be compared to Fig. [3] 
in Sec. 14. 21 where leptogenesis is not required. We see that due to imposing the observed baryon 
asymmetry, the parameter space around J ~ in Fig.[U^ is significantly suppressed and more 
samples are distributed along the upper left-wing (rather than the upper right-wing), while Fig. [9)3 
shows a clear lower bound on ^13. 

Under successful leptogenesis, the correlation between J or Mge with the Dirac CP-phase 5d 
are plotted in Fig.[10l As expected, we see that there are two gaps around 6d = ^ and (Jd = vr , 
respectively. The first gap around do = is more significant. This is consistent with what we have 
observed from Eq. (j5.125p and Fig. [6] earlier, as discussed below Eq. ()5.125p . 



^^S.-F. Ge, H.-J. He, F.-R. Yin, work in preparation. 
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Figure 8: The correlation between 613 and 623 — 45°: all the inputs are the same as Fig. [21 except 
requiring successful leptogenesis in the present figure, for 2000 samples. 
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Figure 9: The correlations of ^13 with Jarlskog invariant J [plot-(a)] and with neutrinoless double 
beta decay observable M^e [plot-(b)], where all the inputs are the same as Fig.[3l except requiring 
successful leptogenesis in the present figure, for 1500 samples. 



Finally, reversing Eq. ()5.123p we can express the leptogenesis scale Mi in terms of baryon asym- 
metry r]Q and other physical parameters, 

47ru2 128(4r2-4rCcos<5D+C^) 



Ml 



(5.129) 



dKjm^ 3[4y-(4r2-4rCcos(5D+C^)^]2 (4r cos (5d - C) sin (Jd ' 

Since the low energy parameters J and Mge in Eqs. (|4.10ip -( |4.102p are also predicted by our soft 
breaking model as functions of (C, ^d) and 5x{= ^13), we see that they must correlate with the 
leptogenesis scale Mi. This is plotted in Fig.[TTl and it shows a robust lower bound on Mi, which 
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is consistent with ()5.127p . as we inferred from Fig. [71 




Figure 10: Under successful leptogenesis, the CP-violation Jarlskog invariant J and the neutrinoless 
double beta decay observable Mee are shown as functions of the Dirac CP-phase 6d , where all 
experimental inputs are varied within their 90%C.L. ranges, for 1500 samples. 
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Figure 11: Correlation of the leptogenesis scale Mi with the low energy Jarlskog invariant J [plot-(a)] 
and the neutrinoless double beta decay observable Mee [plot-(b)] within 90%C.L. ranges, for 3000 
samples. 
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6. Beyond fi — r Symmetry: Dictating Solar Mixing Angle 



In this section we present a general analysis about the derivation of solar mixing angle and about 
how to dictate it from a new minimal hidden symmetry in the seesaw Lagrangian. We have seen 
in Eq. ()4.86p that the solar mixing angle Os (= 612) is fully determined by the elements of Dirac 
mass-matrix m^-, based on the expansion of (r, (^) up to NLO. In Sec. l6.1l we show how the formula 
(j4.86p can be derived in the /i — t symmetric limit. This feature can be proven without expanding 
the seesaw mass-matrix for light neutrinos (Sec. l6.2( ). Actually, we note that 6s is controlled by an 
extra new Z2 symmetry acting on the Dirac mass-matrix shown in Sec. l6.3l 



6.1. Solar Mixing Angle from fi — T Symmetric Mass Matrix 



Here we first derive the solar mixing angle Og as given by Eq. (|2.18bp earlier in the /i — r symmetric 
limit. Let us rewrite Eq. ()2.16p as follows, 



1 



^2aV 



(2-r)rM22 



a{b + c)r a{b + c)r 

- c)2 + 2bcr {I? + c^)r - {b - c)^ 
(6 - c)2 2bcr 



^^0 -Bo ^0 
Co Do 
Co J 



(6.130) 



with an overline on {Aq, Bq, • • • ) to distinguish from the LO elements in Sec. 12. 21 (where further 
expansion in r was made). Thus we have, 

,2 



^0 = 
Co + Do 



2a^ 



2-r)M22 ' 

(6 + c)2 



Co -Do 



a{b + c) 
2-r)M22 ' 
{b - cf 



(6.131) 



(2-r)M22 ' " rM22 

As shown in Sec. 13. H any /i— r symmetric mass matrix can be diagonalized by a two-step rotation, 
because of ^13 = . First, we make a f/23 rotation to partially diagonalize M^^\ with ^23 = 45° (as 
required by the /U — r symmetry), 

/ Ao V2B0 \ 



= C/23(45°)^Mi°^C/23(45°) 



V2B0 Co + Do 
V Co- Do J 

r(0)' 



(6.132) 



where U2z{G) is introduced in Eq. (j3.52cp . Note that the above matrix contains only a non- 

diagonal 2x2 sub-matrix involving {12}-mixing. Second, we further diagonalize M^^^' via the 
rotation U\2 with mixing angle 9s {= ^12) (which is independent of the — r symmetry as noted at 
the end of Sec. 13. II). 



U12 U12 = diag(mi, 7712, mg) , 



(6.133) 
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where U12 is introduced in (|3.52bp . Then, it is straightforward to derivJ^ the mixing angle 6 

2V2B0 



tan 29s 



(6.134) 



as weh as the mass-eigenvalues 
:4o + Co+^o 



m 



1,2 



It 



i[Ao{Co+Do) - 2Bq] 



{Ao + Co+DoY 



rrin 



• 



(6.135) 



In this section, ah mass-eigenvalues are analyzed at the seesaw scale though we have suppressed the 
"hat" notation introduced around the end of SecO 

Note that the minimal neutrino seesaw generally has vanishing determinant and thus. 







detM 



(0) 



^0 {Co + Do) - 2bI {Co - Do) , 



which, for Co 7^ Dq [as indicated in (|6.131|) ]. imposes a condition, 

:4o(Co + :Do) = 2bI. 



(6.136) 



(6.137) 



It is easy to verify this condition by explicitly inputting the LO predictions of (j6.131|) . With (|6.137p 
we can readily deduce the mixing angle from (j6.133p or ()6.134p . 

^0 



V2B0 

as well as the mass-eigenvalues from (|6.134p . 

0, mo = Ao + Co + Do, 



mi 



tan6l., = (6.138) 

= Co-Dq. (6.139) 

We may also derive the above mass-spectrum by a different way. From the trace of M^^^^ or M^^^' , 
we know the sum of three mass-eigenvalues, m-^ + + = ^0 + 2Co , while Eq. ()6.132p gives 
771,3 — Cq — Do ■ So, we have 777^^ + mg = + Co + -Do • The minimal seesaw must have vanishing 
determinant of M^^^ and thus ■m-^^m2m-^ = . Hence, we have, (m^, mg) = (0, ^40 + Co + -Do) or 
{nil, "^2) ~ (^0 + Co + -Do, 0) . But to determine whether 777^ or m2 equals zero will require explicit 
derivation of the rotation matrix Ui2{0s) and eigenvectors, as we did above, which predicts m-^ = , 
corresponding to the NH mass-spectrum for light neutrinos. 

Using the explicit form of the fi—T and CP symmetric mass- matrix ()6.130p -( )6.13ip . we can readily 
deduce the solar mixing angle from (|6.138p . 

V2a 



tan 9s 



b + c' 



(6.140) 



where by a proper phase shift a^^ — )• ai + mr , we can simply convert ()6.140p to ()4.86p or (j2.18bp 
which always ensures tan 9s ^ . The formula (j6.14Up or ()2.18bp also shows that solar angle 9s 
derived from the /U— r symmetric mass-matrix M^^^ is solely determined by the elements of the Dirac 



^The expressions ()6.134|) - H6.135[) agree with tiiat given in Ref. [43] earlier for general /i — r symmetric mass matrix. 
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mass-matrix and is independent of the Major ana mass-matrix M/j . Besides, since M^^'^ is real, 
all CP phases are irrelevant here. Finally, substituting ()6.13ip into ()6.139|1 . we further derive the 
mass-eigenvalues for light neutrinos in the ^ — t symmetric limit, 

2a2 + (6+c)2 (6-c)2 

rm = 0, mo = ; mo = \^ /, , (6.141) 

^ ^ M22(2-r) 3 M22|r| ' ^ ^ 

where we can always ensure > by absorbing the sign of r into the Major ana pahse e^'^^ in U' 
[cf. Eqs. (|3.5U|) and (|3.52ap ]. The mass formula (|6.141|) proves what we gave earlier in Eq. (|2.18ap of 
Sec. 12. II under the /i — r symmetric limit (with the sign choice p = + )■ 

6.2. Solar Mixing Angle Not Affected by Soft Breaking 

In this subsection, we present a more general proof that the solar mixing angle Og is determined by 
the zeroth order fi—T symmetric mass-matrix mI^^ only, independent of any fi—T and CP breaking 
parameter in either SM^ or [cf. (|6.142ap below]^. According to the decomposition method 

in Eqs. (j2.3ip - (j2.32p . we can uniquely split any mass-matrix for light neutrinos into their /U — r 
symmetric and anti-symmetric parts, 

= Mt + M^, (6.142a) 
Mt = mI)^^+6M^, = 5M^, (6.142b) 

where the superscripts ^ and denote the /i— r symmetric and anti-symmetric parts, respectively. The 
/i— r symmetric is due to small 0{y) correction (related to solar mass-squared difference Amg^), 
and 6M^ are induced by the /i— r breaking effects characterized by the two small parameters 5x and 
6a ■ Since neutrino data require (y, 5x, 5a) to be small, it is enough to expand its mass-eigenvalues 
(in Dy) and diagonalization matrix V up to the linear order of (y, 5x^ Sa) for all practical purposes. 
Thus, we write, 

= Dl + 5Dt = (D^^^ + 5Dt) + 5Di , (6.143a) 
V = Vs + 5Va = {Vo + 6V,) + 6Va , (6.143b) 

where Vg corresponds to the ^— r symmetric solution (|3.59p at the end of Sec. 13. H 

Vs = V {Os, 9a = 45°, Ox = 0°|a2 = a^) , (6.143c) 

and 6Va oc {6x, 6a, — ^3) . From (|6.143p - (|6.143cp . we derive the following, 

= (v: + 6V:)D,{V} + 6V^) 

= V:d,V,^ + (v:d,6V^ + 6V:d,V,^ + 6V:DJvA . (6.144) 



previous study [33] analyzed a perturbative violation of a general low energy /i — t symmetric mass-matrix of 
light neutrinos, it also noted that 612 (and mass-eigenvalues) do not receive any /i — r breaking correction. We thank 
an anonymous referee for kindly bringing [44| to our attention. 
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Since the decomposition ()6.142ap is unique, so we can deduce, 



5M"; 



(6.145a) 



(6.145b) 



where 5j denotes the generic NLO parameters under consideration (such as 5x, Sa and y, etc). 

A key point here is to note that the fJ.—T symmetric solution (|3.59|) does not restrict Og and thus 
there is no need to expand 6^ in Vs (despite the ^—t breaking parameters are small). Hence, we can 
determine Og by simply solving Vg from the diagonalization equation ()6.145a|) . In Sec. 16. 11 we have 
used Vg to diagonalize the /i— r symmetric mass-matrix M 



(0) 



M, 



u\C=0 



and found that the formula 



of 9g is determined by m^, only, but independent of M/j. Now with full M*|^^o ; we will prove that 
the formula of Og remains unchanged. 



Inspecting the general form of , 

-- 

we can diagonalize it by the rotation [/23(45°) 
M^' = t/23(45°)^M,^C/23(45°) 



'A Bg 
Cg 

/ 

V 



(6.146) 



A 
V2Bg 




V2Bg 
Cg + D 







Cg-D ) 



(6.147) 



Then the remaining {12} sub-block can be readily diagonalized by the rotation Ui2{ 

A 



with 



tan 6^ 



V2Bg 



(6.148) 



which extends Eq. ()6.138p with the C 7^ effects included. The elements of ()6.146p can be directly 
inferred from the complete seesaw mass-matrix (|2.38ap . without any expansion, 



A 



Bg 



Cg + D 



a\2-X) 


{2-r 


■)(2 


-X) 


{2-r-X)M22 


2(2- 


-r- 


-X) 


a{h+c){2-X) 


(2-r 


)(2 


-X) 


2{2-r-X)M22 


2(2- 


-r- 


-X) 


{b+cf{2-X) 


(2-r 


)(2 


-X) 


2{2-r-X)M22 


2(2- 


-r- 


-X) 



An 



Bo 



(Co + Do) , 



(6.149a) 
(6.149b) 
(6.149c) 



where after the second equality of each equation we have used the ^—t symmetric form in ()6.130p - 
(|6.13ip . The relations (|6.149ap - (|6.149cp prove that the elements {A, Bg, Cg + D) differ from the LO 
results (^0) Bq, Cq + Dq) only by a common overall factor which is irrelevant to the mixing angle 
6g . With these, we explicitly deduce from (|6.148p . 

^0 V2a 



tan 9^ 



V2Bo 



b+c 



(6.150) 



45 



in perfect agreement with Eqs. ()6.138p and ()6.140p . which we derived ear her under the fi—r symmetric 
hmit in Sec. 16. II As noted below (|6.14Up . making a proper phase-shift — )• ai + nvr , we always 
ensure tan > and thus 6s G [0, ^] . Thus we may rewrite (j6.150p as 

where we have used the definition k = (b + c)/a as in (|4.77p . 

This completes our proof that the solar mixing angle 9s is not affected by the /U — r breaking 
and depends on the Dirac mass-matrix only. As a byproduct of the next subsection, we will 
further reveal that this fact actually originates from the elegant feature of a new hidden symmetry 
underlying the seesaw Lagrangian. 

6.3. Beyond /^ — r Symmetry — A New Hidden Symmetry for 6.^^ 

In Sec. l6.Hf6.2l we have revealed that the solar mixing angle Os (= ^12) is fully determined by the 
structure of the Dirac mass-matrix only, independent of the Majorana mass-matrix Mr . In 
particular, Eq. (|6.15ip shows that 9s just depends on a ratio, 

k = (6.152) 

a 

Contemplating this striking feature, we wonder: Is this unique ratio k connected to any new hidden 
symmetry underlying the structure of ? 

The simplest such hidden symmetry we could possibly imagine is a new discrete Z2 symmetry 
which we will denote as Z| to indicate its connection to the solar angle 9s ■ Under this new Z| 
symmetry, we expect the three light neutrinos = (z/g, u^, u^Y^ transform in its 3-dimensional 
representation and the right-handed neutrinos are just singlets (as the above ratio k is independent 
of Mr), 

Ul^TsUl, N^N. (6.153) 

If the hidden is a fundamental flavor symmetry of our theory, it must keep the seesaw Lagrangian 
(j2.2p invariant under the transformation (j6.153p . This means, 

Tlmo = (6.154) 

where the Dirac mass-matrix respects the ^jl—t and CP symmetries in our construction. Thus the 
3x3 unitary transformation matrix Ts should be real and orthogonal, tJt^ = TJ^Ts = X3 , with 
X3 the 3x3 unit matrix. Since € Z| , we must also have T^^ = X3 . Imposing these requirements 
we can solve out Ts from the invariance condition ()6.154p . 

/ 2-/c2 2k 2k 

2k k^ -2 I , (6.155) 
V 2k -2 k"^ 



Tsik) ^ 



2+A;2 
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where k = — ^ as defined in ()6.152p . It is straightforward to exphcitly verify that Ts{k) indeed 
obeys the invariance equation (|6.154|) as weh as the desired conditions, 



2:3, 



(6.156) 



Hence, Ts{k) does form a 3-dimensional real orthogonal representation of the hidden symmetry Z| 
in our seesaw Lagrangian ()2.2p . and the ratio k = is just the group-parameter of . 

We further note that the new Z| symmetry is also generally respected by the low energy seesaw 
mass-matrix ()2.3p . 



T^mj,M^\T^mj,) 



1 rp 



(6.157) 



Revealing this hidden symmetry we can now fully understand why solar mixing angle 6s is not 
affected by the Majorana mass-matrix Mr and the soft fi — r and CP breakings therein. This is 
because, as determined in ()6.15ip . the Og is generally protected by the hidden symmetry via its 
group parameter k , independent of the form of Mn . 

This is an intriguing and essential feature of our soft breaking seesaw model, which poses two 
discrete flavor symmetries C?) Z| in addition to the CP-invariance, among which 1^2^ and CP 
receive small soft breakings from a common origin in Af/j . The exact hidden symmetry Z| dictates 
the solar angle 6s (= ^12)) while the common soft breaking of fi—r and CP connects 6x (= ^13) and 

(= ^23 — 45°) . Under the symmetry group Zi^"^ (X" Z^ and using Eq. (j6.15ip . we can reduce the 
neutrino mixing matrix UiOs, 9a, 9x) in (|3.52bp to 

|fc| 



U{k) = C/(0„45°,O°) 
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Cs 


-Ss 
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V2 
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-V2 

M 



\ 



V2+fe2 ^2{2+fc2) 

1 \k\ _ 



-1 
1 



(6.158) 



We note that this Z| symmetry has a simple geometric interpretation. In the invariance equation 
(|6.154p the Dirac mass-matrix consists of two columns which can be viewed as two vectors, i.e., 



i]^, 1x2; with 



(a, b, c)^ and 



Un 



a, c, 6)"^ in a 3-dimensional coordinate frame. Then, 



the invariance equation ()6.154p just means the eigenvalue equations. 



TsU^ 



u 



1 ' 



TsU2 



(6.159) 



The vectors and U2 must be linearly independent for any realistic , and thus determine a 
plane S. The equation for the plane S reads. 



1 , 

X- -^{y + z) = 0, 



(6.160) 



which is characterized by a single parameter k as defined in ()6.152p . For any other given vector u' , 
so long as it lives in the plane S, it can be expressed as a linear combination of (ui, u^), and thus 
is also an eigenvector of Ts , leading to Tsu' = u' . If u' lives outside the plane S, i.e., disobeys 
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Eq. ()6.160p . then the transformation Tg reflects the vector u' respect to the plane S, as Tgu' . Since 



the plane S has the normal vector n_|_ = {k,—l,—l)/\/2 + k'^ perpendicular to itself, it is clear 
that Tjjn_|_ = —nj_ holds. Then, we can project the vector u' along the directions parallel and 
perpendicular to the plane S, = Uy + u'j_ , so we have 

TgU = Tsiijj + Tsii'_i_ = uj| — , TgU = u, (6.161) 

where uj| lies in the plane S and is thus a linear combination of vectors (u^, Ug) , while u'j_ (x nj_ . 
Hence, Eq. ()6.16ip proves that the operation Tg does reflect the arbitrary vector u' respect to the 
plane S, and u' goes back to itself after applying twice. It is clear that such transformations 
form the 3-dimensional representation of a discrete Z2 group. 

For the case of general seesaw with three right-handed neutrinos M = {Ni, N2, N^)'^ , the Dirac 
mass- matrix nij^ is extended to a 3x3 matrix which, under the /x — r symmetry, takes the form, 

/a' a a \ 

m'jy = b' b c \ (6.162) 
^b' c b I 

where all elements are real due to CP-conservation. Now, to hold the invariance equation (j6.154p we 
just need the flrst column = (a', 6', b')"^ to live in the S plane. Thus the plane equation (j6.160p 
requires, 

4 = = 7^ (6-163) 

2b' b + c k ^ ' 

under which the symmetry Z| holds for general three-neutrino-seesaw, and the solar mixing angle is 

still determined by the group parameter k via Eq. ()6.15ip . We note that the in ()6.162p contains 

two more parameters (a', b') than the minimal seesaw form of , but the Z2 symmetry requires 

only one of them be independent due to the condition (j6.163p . Furthermore, from Eq. (j6.163p we can 

express U3 as a combination of (u^, U2), 

b' 

U3 = -T-, — {u^ + U2), (6.164) 
b + c 

which forces m'^ to be rank-2 and thus det m'j-, = . So, the corresponding seesaw mass-matrix for 
light neutrinos. 

Ml = m'j^M'j^^m'i^ , (6.165) 

must have vanishing determinant because detM^ = (det m^)^(det M^)~^ = for any 3x3 heavy 
Major ana mass-matrix M'^ . It is clear that this feature depends only on . 

Hence, we can state a general theorem: for the Z| symmetric seesaw Lagrangian, the Dirac 
mass-matrix m'j^ is rank-2 at most, enforcing the seesaw mass-matrix M[, to be rank-2, having zero 
determinant and thus one zero-mass-eigenvalu^^. This means that M[, shares a similar feature with 
the minimal seesaw (including two right-handed heavy neutrinos), and they both have, mim,2m^ = . 



^^The possibility of m'o or being rank-1 is already ruled out by the oscillation data where both Aa and As are 
nonzero (Table-[l]), requiring at least two non- vanishing mass-eigenvalues for light neutrinos. 
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Finally, we note that our Z| symmetry allows its group-parameter k to take different values, and 
thus gives different values of the solar mixing angle 6s ■ For instance, to predict the conventional 
tri-bimaximal mixing [35] only needs to assign |A;| = 2 in our construction, leading to tan^^ = 
(6g ~ 35.3°); while we can make another equally simple assignment of \k\ = to generate 
tan 6*8 = I {Os — 33.7°) , which is well within the la range of Og and agrees even better to its central 
value in Table-[TJ 

The hidden symmetry Z| points to a very encouraging direction although it is not yet powerful 
enough to fix the group-parameter k on its own. We expect that at a higher scale the two elegant 
discrete symmetries Zi^'^ (8) can be unified into a larger flavor symmetry which is more restrictive 
and thus fully fixes k. For instance, let us consider a larger group 54 |l6j, the permutation group 
of four objects, which has five irreducible representations {1, 1', 2, 3, 3'} , and assign the left- 
handed neutrinos to 3 . With boldface superscript denoting a given irreducible representation, the 
3-dimensional representations Gj {j = 2, 3) of 54 are, 



/4 i i\ /I 



2 1 2 _ /-y3' 

3 ~3 oil L^:? — 



\ 2 2 

\ 3 3 



1 |, (6.166) 

\ 1 



and the Gf = G3 is non-independent. With these, we can thus identify, 

Gl = -Gf = T^, (6.167a) 
= Gf' = Ts{k = 2), (6.167b) 

where Tj^ is the transformation matrix we defined in ()2.14p . We note that Ts{k = 2) is just a 
special case of our 3-dimensional representation Ts{k) of the hidden Z| symmetry, fixing k = 2 and 
corresponding to the conventional tri-bimaximal ansatz [35], tan^^ = , via our Eq. ()6.15ip . This 
also supports the previous studies of tri-bimaximal ansatz that leads to the ^4 group [3S] • 

Strikingly, our general construction of the hidden symmetry contains the tri-bimaximal mixing 
as a special case (k = 2) and allows deviations from it, e.g., as noted above, another simple choice, 
k = , leads to tan^g = | with a better agreement to the measured value of 6s (Table-[T|). This 
direction will be further explored in a forthcoming publication [47j . 



7. Conclusions 

The oscillation data have provided compelling evidence for the /u— r symmetry as a good approximate 
symmetry in the neutrino sector. The /i— r symmetry requires vanishing mixing angle ^13 and thus 
the Dirac CP-conservation. So, the /i — r breaking term must be small and also serve as the source 
of Dirac CP-violation. On the theory ground, it is natural and tempting to expect all CP violations 
arising from a common origin, implying that both Dirac and Majorana CP violations vanish in the 
/i — r symmetric limit. In this work, based upon these we have conjectured that both discrete — r 
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and CP symmetries are fundamental symmetries of the seesaw Lagrangian (respected by interaction 
terms) and are only softly broken from a common origin. Such soft breaking has to arise from a 
unique dimension-3 Majorana mass-term of the heavy right-handed singlet neutrinos, as shown in 
Sec. l2.2[ This conjecture can hold for general seesaw with three right-handed neutrinos; we have 
chosen the minimal neutrino seesaw [29] in Sec. 2 and shown that the soft ^ — r and CP breakings 
are further characterized by a single complex parameter in the mass-matrix M/j of right-handed 
neutrinos [cf. (|2.29p ]. (In Sec. 6.3 we have proven that the general three-neutrino-seesaw under the 
fJ. — T symmetry and hidden symmetry also predicts a massless light neutrino, sharing the 
same feature as the minimal seesaw.) 

Prom this conceptually attractive and simple construction, we have predicted the soft ^ — t 
and CP breaking effects at low energies. This gives the quantitative correlations between the two 
apparently small deviations ^23" (= ^a) and 0^^— 0° (= 5x)- For any nonzero ^23 ~ ' 
can place a generic lower limit on the mixing angle ^13 in Eq. (j4.80p of Sec. 14.11 This feature is 
demonstrated in Fig. [2] (Sec. 14. 2( 1 without requiring leptogenesis, where ^13 is shown as a function of 
the deviation ^23" 45° in its 90%C.L. range. Fig. [2] also predicts a new upper hound, 6*13 < 6° , for 
full range of ^23 • Adding the successful leptogenesis in Fig. [8] (Sec. l5.2() further predicts a nontrivial 
lower hound, ^13 > 1° , even for ^23 ~ 45° . Fig. [5] and Fig.[S] quantitatively connects the on-going 
measurements of mixing angle ^23 with the upcoming probe of ^13 at the Double-Chooz, Daya Bay, 
T2K and NOi^A experiments, etc. We note that the current measurements of ^23 already show an 
interesting deviation from the maximal value 45° in its central value (due to the subleading effects 
in the oscillation [3 [5]), but still has sizable errors which are larger than that of 612 by a factor of 
3 ~ 4 and are also comparable to that of ^13 [cf. Table-[T] and Eq. Fig. [2] and Fig. [5] reveal 

that a more precisely measured deviation ^23 — 45° will definitely put stronger lower limit on ^13 . 
Hence, our findings strongly encourage experimental efforts to improve the precision of 623 as much 
as possihle. 

We have further derived the low energy Dirac and Majorana CP-violations from a common soft- 
breaking phase associated with fi — T breaking in the neutrino seesaw. In particular, the low energy 
Dirac phase angle 6d and Majorana phase angle (1)23 are predicted in terms of the original soft 
CP-breaking phase angle uj in Mr, by Eqs. ()4.97ap and ()4.96p . respectively. The correlations of 
the low energy Jarlskog invariant J and neutrinoless double decay observable Mee with the mixing 
angle ^13 are analyzed, as depicted in Fig. [3] (without requiring leptogenesis) and Fig.O (with suc- 
cessful leptogenesis). Then, we studied the origin of matter (cosmological baryon asymmetry) via 
leptogenesis, which leads to a robust lower bound on the leptogenesis scale. Mi > 3x10^'^ GeV, as in 
()5.127p of SecO We further analyzed the interplay between the leptogenesis scale Mi and the low 
energy Jarlskog invariant J as well as the neutrinoless double-beta decay observable M^e , shown 
in Fig.[IBi-b. 

Finally, we proved in Sec. 16. Hf6. 21 that the solar mixing angle 612 is independent of the soft fJ.—T 
breaking and only depends on the structure of the Dirac mass-matrix as shown in Eq. (j6.140p or 
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()6.15ip . In Sec. l6.3| we further revealed a new hidden symmetry Z| that dictates the solar mixing 
angle 612 in terms of its group parameter k [cf. (|6.15ip and (|6.155p ]. The new Z| can further 
hold for general three- neutrino-seesaw under the condition (j6.163p . This restricts the 3x3 Dirac 
mass-matrix m'^ to have only one more independent parameter than the corresponding m^, in the 
minimal seesaw. We further proved that the 3x3 Dirac mass-matrix m'^ and the seesaw mass- 
matrix must be rank-2, predicting a zero mass-eigenvalue for ; this same feature is also 
shared in the minimal two-neutrino-seesaw (Sec. 2). In addition, we revealed that the conventional 
tri-bimaximal mixing (TBM) ansatz [l5] is realized as a special case of the group with k = 2 , 
under which (8> [together with the Gil symmetry for leptons in ()2.4p ] is naturally unified 

into a larger group ^4 (supporting previous 54 studies [l6] ) . We note that the Z| symmetry allows 
deviations from the TBM and a different choice of our group parameter k (such as A; = ) can 
give a better agreement to the current data of 9i2 . Further explorations along this direction will be 
given elsewhere |47| . 

Note Added in Proof: 

After the submission of this paper, a newly updated global analysis of solar, atmospheric, reactor 
and accelerator neutrino data for three- neutrino oscillations appeared j48j . which gives the following 
fitted ranges of the low energy neutrino parameters at la {3a) level, 

Am^i = [7.59 ± 0.20 (laeg)] x 10"^ eV^ , 

f [+2.47 ±0.12 (±0.37)1 X 10-3 eV^ (NH), 
Ami = < ' (7.168) 

[ [-2.36 ±0.07 (±0.36)] X 10-3 eV^ (IH), 

012 = 34.5 ± 1.0° (tli:^ , 623 = 42.9+1^: (+^^2^°) , 013 = 6.8+li: (^ 12.8°) , 

where the abbreviations NH and IH stand for the "normal hierarchy" (mi < m2 < m^) and "inverted 
hierarchy" (mi ~ m2 > m^) of light neutrino mass spectra, respectively. Since these updated values 
change very little from the previous global fit summarized in our Table- 1 [3], we find no visible 
correction to all the numerical analyses presented in Figs. 1-11. 
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Appendices 

A. Mass Diagonalization for Right-handed Majorana Neutrinos 

In this Appendix, we perform the mass-diagonalization for heavy right-handed Majorana neutrinos. 
The Majorana mass matrix Mr and its diagonal form are connected by the 2x2 unitary 
rotation Vr , 

V^MrVr = Dr = I , or, Mr = V^DrV^. (A. 169) 

y M2J 

The transformation matrix Vr can be parametrized as 

where {sj^, c^) = (sin^^, cos^^) . Note that there is no independent Dirac CP phase in the 2x2 
unitary matrix Ur and all possible CP-phases are included in the diagonal phase matrices U'j^ and 
U'^ . For notational convenience, we define, 

~ „T „ fM22 M23\ ( 1 (l-r)e^^ \ 

^UMe-. Y ^^^^^^^ 

Then the diagonalization equation (|A.169p takes the following form, 

UIMrUr = Dr, or, Mr = U*rDrUI. (A.172) 

Using Eq. (|A.172p we can reconstruct Mr , 

4Mi + 4M2 CRSR (Mi - M2 , , 
Mr = \ ""^^ ^ ' (A.173) 

, CRSR Ml - M2 4Mi + c\M2 
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Thus, comparing the right-hand-sides (RHS's) of Eqs. (|A.171ap and (|A.173p we can generaUy solve 
the masses Mi and M2 , 



Ml = - ( M22 + M33] + 

M23 



M2 = J f M22 + M33 

2 V / ZCnS 



(A.174a) 
(A. 174b) 



as well as the rotation angle 0^ . 

tan 28 f> = 



2M; 



23 



2(1 -r) 



(A.175) 



As tan 20^ is real, the imaginary part on the RHS of Eq. ()A.175p must vanish. This allows us to 
resolve j3 and 9^ from ()A.175p . 



/3 = -Csina; + 0(a 
?H = ±---Ccosw. 



(A.176a) 
(A.176b) 



The sign of it^ will be fixed from Mi < M2 , as required by ()2.23ap . Substituting ()A.176bp back 
into (|A.174p . we deduce. 



M- 



M22e*('^i+^) 
M2 = M22e^(2^2+^) 



1 T 1 ± r - |e^^ + 0{C^] 



(A.177a) 
(A. 177b) 



We see that requiring Mi ^ M2 as in (j2.23ap will peak up the minus sign for j on the RHS of 
Eq. (|A.176bp . Thus we further simphfy ()A.177aP - (|A.177bp as 

1 



Ml ~ M22e^(2Ti+/3) 

M2 ~ M22e^(2T2+/3) 



2 - r - -Ce^'^ 
2^ 



(A. 178a) 
(A. 178b) 



The mass eigenvalues Mi and M2 are real and positive, so taking the absolute value on both sides 
of (|A.177p we can derive Mi and M2 under the expansion of r and The formulas turn out to fully 
coincide with that in Eq. ()4.106p which we derived earlier. Finally, requiring the RHS of Eq. ()A.178p 
be real and noting f3 = 0{() , we deduce the phase angles, 

^( sinuj 



sin 27]^ ~ 



72 



1 



[r^ —r( cos a; -|- 1^2] 2 
Csino; = 0(0- 



cos 27^ ~ 



r — ^( cos uj 



[r 2 — r C cos w -|- 4 



(A.179a) 
(A.179b) 



As expected, here we see that the three phase angles (71, 72, /3) would vanish as Csincj — )• . 
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B. Positivity Constraints from Baryon Asymmetry 77^ 

In this Appendix we systematically derive the positivity constraints from the baryon asymmetry 
{r]^ > 0) , on the CP-phase 5d of our model, as given in Eqs. ()5.125p of Sec. l5.2l For convenience, 
we rewrite Eq. ()4.99p as 



c 



(B.180) 



where we have used r± to denote the two sign-combinations in the brackets. With the above formula 
we further derive the positivity condition ()5.124|) of i]^ as. 



Csin^D ^cos25D±2cos(^Z)y^^y -sin^^D j > 0, (B.181) 

which reflects the requirement of the baryon asymmetry 7]^ > . Below we will solve this inequality 
and derive the physical regions for the Dirac CP-phase 6d as allowed by the successful leptogenesis. 

B.l. Solution of Eq. ( ]B.181[) for r = = r_ Branch 

The case of r = r_|_ = r_ corresponds Csin5/) cos 25/) > 0. Since our convention always holds 
C > , the condition ()B.18ip reduces to 

sin(5£)Cos2(5z) > 0, (B.182) 

which requires 

, („, ^) u (5^, ,) U II) . (B.183, 



B.2. Solution of Eq. (IB.lSlh for r = ^ r_ Branch 

The case of r = r , corresponds to -|- sign in (IB. 1811) . WithC > 0, the condition (IB. 1811) reduces to 



sin^D ^cos2(^Z) -h2cosfey||^^ -sin^Ji) j >0. (B.184) 

The defined range of 6d is [0, 2tt) , but the inequality ()B.18ip requires sin 5/) ^ and thus excludes 
6d = 0, ■ So, we will analyze the ranges 6d G (0, vr) and 6d G (vr, 2tt) , respectively. 

B.2.1. For the Range So G (0, vr) 

For 6d G (0, vr) , the condition (|B.184p b ecomes, 



cos26d > -2cos(5dW-^^ -sin^Js , (B.185) 
V 16 



which can be resolved as below. 
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For the sub-range (5d G (O, ^) , we have cos2(5z) > and —cos6d < 0, which always hold 
the condition (|B.185p . 

For the sub-range 6d ^ {j, ^) , we have cos 26d < and — cos6d < , leading to the solution 
from (1B.185P . 

cos^^^ > 4479- (B.186) 

For the sub-range So G ^) , we have cos 25/) < and — cos^d > 0, showing that the 
condition (|B.185p has no solution. 

For the sub-range 5d G ' have cos 25/) > and —cos6d > 0, leading to the 

solution from ()B.185|) . 

4 



4 S 

cos^ 6d < (B.187) 

st y C 



B.2.2. For the Range So G (tt, 27r) 

For 5z) G (vr, 27r) , the condition ()B.184h becomes, 



/ „4 y2/'2 

cos2(5d < -2cos(5dW— ^-r^ - sin^J/) , (B.188) 
V 16 6^ 

which can be resolved as below. 

• For the sub-range 5/) G (vr, we have cos 25/) > and —cos (5/) > 0. From the condition 
()B.188p . we deduce the solution just as in (|B.186p . 

• For the sub-range G we have cos 25/) < and — cos 5/) > 0, which always hold 
the condition (|B.188I) . 

• For the sub-range 5/) G ^) , we have cos 25/) < and — cos 5/) < . From the condition 
()B.188p . we deduce the solution just as in (|B.187p . 

• For the sub-range 5/) G 27r), we have cos 25/) > and —cos 5/) < 0, showing that the 
condition (|B.188P has no solution. 



B.3. Solution of Eq. (IB. 1811) for r = r 7^ r+ Branch 

The case of r = r_ corresponds to — sign in (|B.18ip . With > in our convention, we reduce the 
condition (|B.18ip to 



sin5/) ^cos25/) -2cos5/)y^^^ -sin2 5/) j >0, (B.189) 
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which excludes = 0, vr. So, we wiU analyze the ranges do G (0, vr) and do G (tt, 27r) , respec- 
tively. 

B.3.1. For the Range Sd G (0, vr) 

For 6d G (0, vr) , the condition (IB.189P becomes, 



J 4 2/-2 
— -sin^f^D , (B.190) 



16 6t 



which can be resolved as below. 



• For the sub-range 6d G (O, j), we have cos2(^£i > and cos6d > 0. From the condition 
(IB.190p . we deduce the solution just as in ()B.187|) . 

• For the sub-range 6d £ (f , f), we have cos26d < and cos 5/) > 0, showing the condition 
(|B.190P has no solution. 

• For the sub-range So G we have cos25£) < and cosSd < 0. From the condition 
(IB.190p . we deduce the solution just as in ()B.186|) . 

• For the sub-range 6b ^ (x' '''")' have cos 25/) > and cos (5/5 < 0, which always hold the 
condition (|B.190p . 

B.3.2. For the Range So G (vr, 27r) 

For 6d G (ti", 27r) , the condition ()B.189p becomes, 



/ „4 y2/-2 

COS 25c < 2cos5dW-^^ -sin^^i? , (B.191) 
V 16 54 

which can be resolved as below. 

• For the sub-range 6d S (tt, we have cos25/p > and cos 5/) < 0, showing that the 
condition (|B.19ip has no solution. 

• For the sub-range 5d G (x' T")' have cos 2(5/) < and cos5d < 0. From the condition 
(|B.19ip . we deduce the solution just as in (|B.187p . 

• For the sub-range 6d G x)' have cos 25/) < and cos5d > 0, which always hold 
the condition (IB.19ip . 

• For the sub-range 5d G (x' ^^)' have cos 25/) > and cos 5/) > 0. From the condition 
(|B.19ip . we deduce the solution just as in (|B.186p . 



56 



Finally, we can summarize all solutions to the positivity condition (jB.lSOp into a compact form, 
for r = r I = r_ , r = r , 7^ r„ and r = r_ 7^ r, , 



for r 



.„.,0,3u(|„)u(|,f); 



for r = r, ^ r- , 



IT 



So € ^0, - 

for r = r_ 7^ r+ 



vrvrx /Stt \ / SvrN 



T' T 



/Svr 77r\ 



Svr \ / Svr Svr , 
— , vr U — , — U 
4 ' / V 4 ' 2 



Svr 7vr 

T' T 



U ( ^, 2vr 



(B.192a) 



(B.192b) 



(B.192c) 



where the regions with subscript ">" require imposing the condition (|B.186p and those with subscript 
"<" are constrained by the condition ()B.187p . The above constraints just confirm what we have 
presented in Eq. (|5.125p of Sec. 15. 21 
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